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Ab-initio density functional calculations on explicitly doped La2-iSra;Cu04 find doping creates 
localized holes in out-of-plane orbitals. A model for cuprate superconductivity is developed based 
on the assumption that doping leads to the formation of holes on a four-site Cu plaquette composed 
of the out-of-plane Ai orbitals apical O Pz, planar Cu d322_^2, and planar O This is in contrast 
to the assumption of hole doping into planar Cu d^2_y2 and O pa orbitals as in the t-J model. 
Allowing these holes to interact with the d^ spin background leads to chiral polarons with either 
a clockwise or anti-clockwise charge current. When the polaron plaquettes percolate through the 
crystal at a; ~ 0.05 for La2-iSra;Cu04, a Cu d^2„y2 and planar O pcr band is formed. The computed 
percolation doping of a; ~ 0.05 equals the observed transition to the "metallic" and superconducting 
phase for La2-2,SriCu04. Spin exchange Coulomb repulsion with chiral polarons leads to D-wave 
superconducting pairing. The equivalent of the Debye energy in phonon superconductivity is the 
maximum energy separation between a chiral polaron and its time-reversed partner. This energy 
separation is on the order of the antiferromagnetic spin coupling energy, Jdd ~ 0.1 eV, suggesting a 
higher critical temperature. An additive skew-scattering contribution to the Hall effect is induced 
by chiral polarons and leads to a temperature dependent Hall effect that fits the measured values 
for La2-a;Sra;Cu04. The integrated imaginary susceptibility, observed by neutron spin scattering, 
satisfies lu/T scaling due to chirality and spin-flip scattering of polarons along with a uniform dis- 
tribution of polaron energy splittings. The derived functional form is compatible with experiments. 
The static spin structure factor for chiral spin coupling of the polarons to the undoped antiferro- 
magnetic Cu dP spins is computed for classical spins on large two dimensional lattices and is found 
to be incommensurate with a separation distance from {n/a,n/a) given by 5Q ~ {2n/a)x where x is 
the doping. When the perturbed — band energy in mean-field is included, incommensurability 
along the Cu— O bond direction is favored. A resistivity ~ yA'+i arises when the polaron energy 
separation density is of the form ~ A'' due to Coulomb scattering of the — y^ band with polarons. 
A uniform density leads to linear resistivity. The coupling of the x^ — y^ band to the undoped Cu 
d^ spins leads to the angle resolved photoemission pseudogap and its qualitative doping and tem- 
perature dependence. The chiral plaquette polaron leads to an explanation of the evolution of the 
bi-layer splitting in Bi-2212. 

PACS numbers: 71.15.Mb, 71.27.-l-a, 74.25. Jb, 74.72.-h 



I. INTRODUCTION 



It is generally assumed that the relevant orbitals for 
understanding high temperature cuprate superconduc- 
tivity arise from holes on planar Cu d^^-y^ and O Pa 
orbitals. The t-J modeli and its generalization to the 
three band Hubbard modeP are believed to be the cor- 
rect Hamiltonians for understanding these materials. Ex- 
tensive work since the original discovery^ has not led to a 
complete understanding of the properties of the cuprates, 
despite the rich physics contained in such a simple Hamil- 
tonian. 

In this paper, we assume doping creates polarons com- 
posed of apical O pz hybridized with Cu d322_^2 and pla- 
nar O Per that form localized chiral states in the vicinity 
of the dopant (Sr in La2-a:Sra;Cu04, for example). The 
polaron orbital is spread over the 4-site Cu plaquette near 
the Sr and is stabilized in a chiral state due to its interac- 
tion with the antiferromagnetic cP spins on the undoped 
Cu sites. This is similar to prior work^i^i^i^ suggesting 
chiral spins states arise from doping except, in our case, 
the polaron is formed from out-of-plane orbitals. 



As the doping is increased, the chiral polarons even- 
tually percolate through the crystal. We assume a Cu 
d,j.2_y2 and O Pa delocalized band is formed in the per- 
colating swath. This leads to our Hamiltonian of a delo- 
calized Cu d^2_y2 band interacting with chiral plaquette 
polarons and localized d^ antiferromagnetic spins on the 
undoped Cu sites. 

For low dopings, momentum k is not a good quantum 
number because the — /pa band is formed on the 
percolating swath. This leads to broadening observed in 
angle-resolved photoemission (ARPES) measurements. 
As the doping is increased, k becomes a better quantum 
number. 

With increasing doping the 4-site chiral polarons crowd 
together in the crystal and several changes occur. First, 
the apical O and single Cu closest to Sr is doped with 
Pz instead of the four Cu's of the plaquette.^ Second, 
the reduction of undoped d^ spins decreases the energy 
difference between a polaron state and its time-reversed 
partner. Third, the number of — jpa band electrons 
increases. 

In our model, the superconducting d-wave pairing is 
due to the Coulomb spin exchange interaction of the 
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— y^ band with chiral polarons where the Debye energy 
in phonon superconductors is replaced by the maximum 
energy difference of a polaron with its time-reversed part- 
ner (the polaron with flipped chirality and spin). This 
leads to an overdopcd phase where superconductivity is 
suppressed. 

Calculations in this paper of the doping values of 
La2-2;Sr2.Cu04 and YBa2Cu306+5 for percolation of po- 
laron plaquettes and the formation of an — band 
are x « 0.05 and 5 m 0.36. Percolation of single doped 
apical O pz and Cu described above is a: ~ 0.20. 
These numbers are close to known phase transitions in 
La2-a;Sra;Cu0 4^°i^^ {x ~ 0.05 for the spin-glass to super- 
conducting transition and x « 0.20 for the orthorhombic 
to tetragonal transition) and YBa2Cu306+(5^^ {5 ~ 0.35 
for the antiferromagnetic to superconducting transition). 

Chiral polarons couple to the Cu (P spins on the 
undoped sites and distort the antiferromagnetic order 
leading to incommensurate magnetic neutron scattering 
peaks. The charge current of the polaron induces a chiral 
couphng of the form ± Jch[Sz • (S^i x 8^2)}^'^'^ where Sz 
is the polaron spin and the subscripts dl and d2 repre- 
sent Cu x^ — spins at adjacent sites. The sign of the 
interaction is determined by the chirality of the polaron. 
This term is in addition to an antiferromagnetic coupling 
between the polaron spin and a neighboring spin, Jdz, 
and the d^ — d^ spin coupling, Jdd- 

We have performed energy minimizations on large lat- 
tices of classical spins doped with chiral plaquette po- 
larons over a range of coupling parameters to compute 
the static spin structure factor. These calculations are 
similar to previous computations of the correlation length 
and incommensurability due to chiral plaquettes^"* us- 
ing the Grempel algorithmic to search for a global min- 
imum. A neutron incommensurability peak consistently 
appears on a circle in k-space centered at (tt/c, tt/o) with 
radius « (2TT/a)x. This result is missing the kinetic en- 
ergy perturbation of the d^2_y2 band. Computing this 
contribution in mean field selects the incommensurate 
peaks along the Cu-0 bond directions in accord with 
experiments ji0ii5ii6,i7 

If the energy difference between a chiral state and 
its time-reversed partner, where the spin and chiral- 
ity are flipped, is uniformly distributed over an en- 
ergy range larger than the temperature, then the dy- 
namical magnetic response of the polarons satisfies uj /T 
scaling J ^ ° ■ ^ ' ^ ^ ^! ^ ^' ^ ^ Since the polarons are ran- 
domly distributed throughout the crystal with different 
undoped d^ environments, the probability distribution of 
the energy separation of these states may be approxi- 
mately uniform. 

There are four possible orbital state symmetries for a 
polaron delocalized over a four Cu plaquette. They are S, 
Dxy, and Px' ± iPy' ^ where the last two states are chiral. 
x', y' refer to axes along the diagonals. Coulomb scat- 
tering of Cu x"^ — y^ band electrons with polarons leads 
to a linear resistivity in the case of a uniform energy dis- 
tribution of the energies of the four polaron states. This 



may be uniform for the same reasons discussed above 
for neutron scaling. Any non-uniformity of the energy 
distribution spectrum makes the resistivity non-linear. 

Spin exchange Coulomb scattering of an x^ — y^ Cooper 
pair (k ti ~k |) with a chiral polaron P^' ± iPy' and spin 
s into the time-reversed intermediate state P^' T 'i'Py' a-nd 
spin — s leads to an anisotropic repulsion that is peaked 
for scattering of a Cooper pair with k near (±7r, 0) to k' 
near (0,±7r). There are two necessary conditions to ob- 
tain d-wave superconducting pairing. First, time-reversal 
symmetry must be broken such that P^' + iPy' and spin s 
is not degenerate with P^i — iPyi and spin —s. The max- 
imum energy separation of these two polarons replaces 
the Debye energy in phonon superconductivity. Second, 
the polaron must be spread out over more than one site 
so that phase differences in the initial and final x'^ — y^ 
band states can interfere. A single site polaron would 
lead to an isotropic repulsion and no superconductivity. 

In zero magnetic field, there is an equal number of po- 
larons of each chirality. A magnetic field creates more 
polarons of one chirality than the other. An — 
band electron scattering from a chiral polaron is skew- 
scattered2^iS£i^i2^ due to second-order Coulomb repul- 
sion with a polaron where the polaron orbital changes in 
the intermediate state. This leads to an additive skew- 
scattering contribution to the Hall effect proportional to 
the difference of the number of "plus" and "minus" po- 
larons. For high temperatures, the difference is ~ 1/T. 

For the hole-doped cuprates, the polarons are holes. 
The Coulomb matrix element is negative, U < 0, since 
the change in the Hamiltonian amounts to the removal 
of a Coulomb coupling. Although we have not identified 
the nature of the polaron in the electron-doped system 
Nd2-a;Ce2;Cu04, the same argument makes U > 0. 

The skew-scattering contribution is derived and com- 
puted for reasonable values of the parameters. It is found 
the sign change between the hole-doped and electron- 
doped cuprates appears due to the sign change of U. The 
magnitude of the skew-scattering term is shown to be 
large enough to account for the experimental data. The 
derived functional form for the temperature dependence 
is shown to fit data^^ on La2-a;Sra;Cu04. To our knowl- 
edge, the only explanation for this sign difference between 
the hole and electron-doped cuprates arises from the 
additional (7r,7r) nesting of the Nd2-a;Ce2;Cu04 Fermi 
surface.'^'' 

The d^ undoped spins interact with the x'^ — y"^ band 
electrons. They induce a coupling of a state with mo- 
mentum k to k ± Q where Q ~ (tt, tt) is the incom- 
mensurate peak momentum. This leads to an ARPES 
pseudogapi^ii^i^^i^ The strength of the d^ antiferromag- 
netism decreases with increasing temperature making the 
pseudogap close with temperature. At low doping, there 
are more undoped d^ spins and the coupling to the x^—y^ 
band electrons is larger than the coupling at higher dop- 
ing. The pseudogap increases with decreasing doping 
while Tc is reduced. The couplings leading to the pseu- 
dogap and superconductivity are different in our model. 
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The outline of the paper is as follows. In Section [TTl 
the experimental and theoretical arguments for the ex- 
istence of Pz holes with doping are examined. In par- 
ticular, experiments considered to establish the valid- 
ity of the t-J model'^'5^^- and preclude any substantial 
out-of-plane character are addressed.— Section |TTT] de- 
fines the chiral plaquette polarons. Section Hvl calculates 
the percolation phase transitions and compares them to 
the La2-2;Sr2;Cu04 and YBa2Cu306+5 phase diagrams. 
Section |V] describes classical spin calculations of the neu- 
tron structure factor including the effect of the (P and 
polaron spin incommensurability on the kinetic energy 
of the band — electrons. Incommensurate peaks 
along the Cu— O bond direction are obtained. The po- 
laron magnetic susceptibility is calculated, assuming a 
uniform probability distribution of polaron energy level 
separations, and is shown to satisfy uj/T scaling. In 
Section I VII the Coulomb interactions of — band 
states with chiral plaquette polarons is examined to de- 
termine the possible superconducting pairing channels. 
The spin exchange interaction leads to an anisotropic re- 
pulsion of the form sufficient to create a d-wave gap with 
nodes. Section fVIII describes the resistivity and Hall ef- 
fect due to Coulomb interactions with chiral polarons. If 
the distribution of energy separations of polarons states 
with different symmetries is uniform, then the resistiv- 
ity is linear. A magnetic field produces a difference in 
"up" and "down" chiral polarons leading to an additive 
skew-scattering contribution to the ordinary band Hall 
effect with a temperature dependence consistent with 
measurements. The magnitude and temperature depen- 
dence of the skew-scattering is calculated. Section lVIII Al 
describes our model of the ARPES pseudogap and its 
doping and temperature dependence. Section IVniB] dis- 
cusses the doping and temperature dependence of the 
bilayer splitting observed in ARPES on Bi-2212. Sec- 
tion |IX] discusses the NMR data of Takigawa et al'^* that 
is assumed to be strong evidence for a one-component 
theory because of the similar temperature dependencies 
of the Knight shifts of planar Cu and O in underdoped 
YBa2Cu3 06.63- We argue qualitatively that these results 
are compatible with our model. Section |X] presents our 
conclusions. 



II. EXISTENCE OF Ai HOLES 
A. Ab-initio Calculations 

Becke-3-Lee- Yang-Parr (B3LYP) is a three parame- 
ter hybrid density functional that includes 20% exact 
Hartree-Fock exchange i^i^i^i^i^ Its success has ex- 
tended beyond its original domain of parametrization to 
include the thermochemistry of compounds containing 
transitions metals 

Several years agO)^ we performed ab-initio periodic 
band structure computations using the spin unrestricted 
B3LYP functional on undoped La2Cu04 and explic- 



TABLE I: Ab-initio La2Cu04 bandgap results. LSDA stands 
for local spin density approximation, SIC is self-interaction 
correction, HF is Hartree-Fock, and UB3LYP is unrestricted 
spin B3LYP. The last line is the experimental gap. 



Metnoa 


Bandgap (eV) 


Author 


LSDA 


0.0 


Yu et al.— (1987) 


LSDA 


0.0 


Mattheis3*a (1987) 


LSDA 


0.0 


Picket1»SS (1989) 


SIC-LSDA 


1.0 


SvaneS (1992) 


LSDA-hU 


2.3 


Anisimov et al.^^ (1992) 


SIC-LSDA 


2.1 


Temmerman et al.-^ (1993) 


LSDA+U 


1.7 


Czyzyk et al,^ (1994) 


HF 


17.0 


Su et al.-^ (1999) 


UB3LYP 


2.0 


Perry et al.i^ (2001) 


Experiment 


2.0 


Cinder et al.i^ (1988) 



itly doped La2_a;Sr2;Cu04. For the undoped insulator, 
the antiferromagnetic insulator with the experimental 
bandgap of 2.0 eV was obtained4^ 

Prior to this calculation, the insulating state had been 
obtained by extending local spin density (LSD) compu- 
tations, which yielded zero gap or a metal, to approx- 
imately incorporate the self-interaction correction not 
accounted for in this functional. Table |T] summarizes 
chronologically corrections to the initial LDA results and 
their computed bandgaps. 

Our result showed that an off-the-shelf functional with 
an established track record— i^i^i^ for molecular sys- 
tems could reproduce the results of more elaborate LDA 
corrections. 

In addition, wc found the highest occupied states to 
have more out-of-plane orbital character (apical O pz and 
Cu z^) than obtained by LDA. Svane^'^ also made this 
observation in his self-interaction corrected (SIC) com- 
putation. 

In a second paper ^2 we explicitly doped La2Cu04 with 
Sr to form supercells of La2-a:Sra;Cu04 at special dopings 
of X = 0.125, 0.25, and 0.50. We found an additional hole 
was formed for each Sr atom that localized in the vicinity 
of the dopant of apical O Pz, Cu , and an Aig combi- 
nation of planar O Pa character. The Cu sites split into 
undoped and doped sites. The undoped sites had a 
x^ — y^ hole and the doped sites were still predominantly 

with a mixture of x^ —y^ and hole character. There 
was corresponding hole character on the neighboring O 
atoms in and out of the plane with the appropriate Bi^ 
and Aig symmetries. This led us to argue that out-of- 
plane hole orbitals are a generic characteristic of cuprates 
and must be considered in developing theories of these 
materials. 

At the time, B3LYP had an established track record 
with molecular systems, but its use for crystal band 
structures was in its infancy. This is likely due to the 
difficulty of including exact Hartree-Fock exchange into 
periodic band structure codes. 
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FIG. 1: La2Cu04 with one doped Sr atom. Ab-initio 
calculations^ find an Ai hole localized above the Sr with hole 
character on the apical O Pz, Cu z^, and planar O Po- orbitals. 
The pz character above the doped Cu is smaller than the pz 
below the Cu leading to an Ai state rather than Aig. 



Since the appearance of our doped Sr work, it has been 
found that B3LYP does remarkably well at obtaining the 
bandgaps of insulatorsi ^^i^^'^°i^^ Hybrid functionals ap- 
pear to compensate the overestimate of the gap from 
Hartree-Fock with the underestimate arising from local 
density and gradient corrected functionals. Thus, we be- 
lieve density functionals have established the existence of 
non-planar hole character in La2-a;Sr3;Cu04. 

For La2-2:Sr2,Cu04, there are five Cu sites in the vicin- 
ity of a Sr atom in two distinct Cu02 planes. The Sr is 
centered over four Cu in a square plaquette. The fifth 
Cu couples to the Sr through the neighboring apical O 
between them as shown in figure [T] The hole state com- 
posed of apical O Pz , Cu , and planar O Pa as shown 
in figure [T] appeared with Sr doping. 

The polaron state with the hole delocalized over two 
diagonally opposed Cu in the four Cu plaquette is higher 
in energy in our ab-initio calculation by 0.57 eV for each 
Sr or 0.071 eV for each formula unit Lai.875Sro.i25Cu04. 
The value of 0.57 eV is an upper bound since our ge- 
ometry optimizations only allowed the apical O sites to 
relax. The polaron localizes on two Cu sites due to spin 
exchange coupling with the — hole and the anti- 
ferromagnetic spin ordering of the — holes in our 
periodic supercells. In this paper, the hole state in figure 
[1] obtained from our ab-initio calculations is not taken to 
be the correct polaron. Instead, we postulate Sr doping 
leads to chiral polarons over the four plaquette Cu atoms 
shown in figure [51 This is discussed in section IIIII 



This paper explores the consequences of the assump- 
tions that Sr doping causes holes to appear in Cu four- 
site plaquettes and the most stable configurations are the 
chiral states Px' ± iPy' ■ 

From an ab-initio standpoint, our first assumption is 
plausible for La2-a;Sra;Cu04, but unproven. This may 
be due to the limitation of the special periodic super- 
cells that were chosen out of necessity to perform the 
computation, the restrictive geometry relaxation for the 
plaquette polaron, or it may be a limitation of the B3LYP 
functional. 

For YBa2Cu306+5, we do not have ab-initio proof for 
doping of four-site Cu plaquettes in the Cu02 plane ei- 
ther. In fact, any polaron plaquettes would likely be in 
the yz plane where the Cu— O chains are along the y-axis 
and the z-axis is normal to the Cu02 planes. One way 
in which polaron plaquettes can arise is when two adja- 
cent Cu— O chains each have an occupied O separated by 
one lattice spacing along the x-axis (perpendicular to the 
direction of the chains). In other words, the two chain 
O reside in neighboring chains with minimum separation 
between them. This may create four-site polarons on the 
two Cu02 planes above and below the two O atoms. For 
this paper, the chiral plaquette polarons in figure [2] are 
assumed. 

The second assumption, that the polarons are chiral, is 
true for a localized polaron interacting with an infinite cP 
antiferromagnetic lattice in two dimensions (2D}^i^i^iii^ 
by mapping the 2D Heisenberg antiferromagnet to a con- 
tinuum model and analyzing the effective Hamiltonian 
arising from a path integral formulation. These papers 
did not specifically consider an out-of-plane hole, but the 
analysis is applicable in our scenario. This is discussed 
further in section [1111 



B. Experiment 

Resonance circular dichroism photoemission investi- 
gating the spin of the occupied states near the Fermi 
leve P^i"^^ find a preponderance of singlet occupied states 
just below the Fermi level in CuO and Bi2Sr2CaCu208-i-5 
(2212). These results are considered strong evidence in 
favor of the correctness of the t-J model. In particular, 
it is expected that out-of-plane z^, and Ai^ p^ would 
lead to triplet occupied states near due to exchange 
Coulomb coupling to the orthogonal x^—y^ orbital. Since 
the prima facie evidence is against our proposal, we re- 
view the measurement and its interpretation. 

We show our assumption of a delocalized x^ — y^ band 
on the percolating out-of-plane polaron doped Cu sites 
leads to a null effect for resonance absorption on these 
sites. This arises because a delocalized x^ — y^ band elec- 
tron spin has no correlation to the polaron spin. Thus, 
the experiment measures the spin of the highest occupied 
states on the undoped Cu d^ sites where it is expected 
the first holes would be created in Big combinations of 
ligand planar O Pc orbitals that form a singlet with the 




FIG. 2: Orbital schematic of chiral polarons postulated in our model. We assume these chiral hole states are the most stable 
due to interactions with the undoped (f Cu lattice spins. 



x'^ — hole (the Zhang- Rice^ singlet) . 

The idea behind the dichroism experiment is to use 
circularly polarized incident soft x-rays tuned to the Cu 
L3 (2^3/2) white hue energy (« 931.5 eV). The incident 
x-rays induce the photoabsorption transition 2p^3(fi + 
huj 2p^d^'^ that Auger decays to an ARPES final state 
2p^(P -\- e. The spin-orbit energy separation of the core- 
hole 2pi/2 and 2^3/2 states is sufficiently large (~ 20 eV) 
to guarantee the intermediate state is a j = 3/2 core-hole. 

By monitoring the outgoing electron energy and spin 
along the incident photon direction for each photon he- 
licity, (j^ and cr", the total spin of the final Cu d® is 
obtained. 

In the Bi-2212 experiment, ■^^ the photon is incident 
normal to the Cu02 planes. The analysis below is 
for a normally incident photon. The transition rates 
are slightly different for the CuO case?^ where a poly- 
crystalline sample was used. 

Consider a Cu initially in the d^ state \2p^\'idz2 
; 3^2.2 _j,2 ] -l))? where our notation shows the 

occupied electrons. The dxy^ d^z, dyz orbitals are always 
doubly occupied and are omitted in the wavefunction for 
convenience. The Zd^-i orbital is doubly occupied and the 
single — electron is in a spin state along a direction 
that may be different from the incident photon direction. 
It is represented as a linear combination of t and | along 
the incident photon direction with p -I- p = 1. By 
summing over all helicities and exiting electron spin di- 
rections, the photocmission becomes independent of the 
initial direction of the — electron as shown below. 

Writing x, y, and z for the angular part of the Carte- 
sian variables, x/r, y/r, and z/r, the relevant wave func- 
tions and photon polarization operators may be written 



as c?322_r2 = C2(3z^ — 1), dx2-y2 = \/iC2{x^ — y^), and 
Yi±i = Ci(l/V^)(i ± y), where Ci = ^Jfl^ and C2 = 
y'S/l&TT. The mod-squared matrix elements for reso- 
nance absorption of \2p^]idz2 Ti; 3(ia;2_j^2 (Af f A-Ai [)) 
to the intermediate 2^3/2 core-hole states are, 

|2/;34^ Ti;34^_,2(AT T+Ax [)) 




1 (Ch 

6 1^2 



(1) 



|2/;3d,2 U■Mx--v<A^^ +Ai [)) 



|p^:f,-|)di" 



n5 . 3 



^ IP- ;,,-!)<;■» = L\A,f J (i Vft 

where Tito'^ are positively and negatively circularly polar- 
ized photons. 

The Auger scattering rates of the four 2^3/2 interme- 
diate states where one — electron fills the 2^3/2 core 
hole and the other is ejected are. 



2'2/ 



\2p')\dz2]i)+e]=\V\\ (3) 



2'2/ 



|2p6)|4. n)+eT=(|)|l^P, 
|2/)|42 n)+ei= {r)\V\\ 



(4) 



-jio 



|2p6)|4. U)+eT=(i)|Vp, 
|2p6)|42n)+ei=(|)|y|2, 
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(6) 



where \V\'^ is the Auger matrix element. The is dou- 
bly occupied making the state a singlet. There are 
analogous matrix elements if the Auger process scatters 
the two electrons instead of — y^and also if the final 
d^ is composed of one electron in and one in — 
in a singlet configuration. 

The total scattering rate is given by the products 
through the various intermediate states. Using the 
conventio n^^i^^ cr"*" 1, a~ 1, J,, and a~ | to repre- 
sent a positively circularly polarized photon ejecting an 
electron with | spin etc, the scattering leaving a singlet 

final state is, 



lAr + gi^irli^r, 



a- i 



\Ai\']\V\' 



The total parallel and anti-parallel scattering is. 



TT= (a+ T +a- i) 
ti= {a+ i +a- T) 
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\V\', 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



where we have neglected the (C1/6C2) from equations [T] 
and [2] since it cancels out when we evaluate the polariza- 
tion defined below. These two sums are independent of 
the starting spin orientation of the — electron. The 
"polarization," defined as a ratio (tt — Ti)/(TT + Ti) = 
5/6 for pure singlet d^ states. 

There are three possible triplet d^ spin states. There 
is one electron in x'^ — y^ and z^. The scattering from the 
intermediate d^'~' state with a 2^3/2 core- hole to triplet d^ 
is given by. 



3 1 
2' 2 
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|2p6)|Tmi; 



ei-2|T/|2, 
+ eT- \V\'', 



(13) 



|2/)|U>+ei= (I) -21^^12, 
|2p6)|tt)+et=(i) -21^^12, 



(14) 



'---)d^ 

2' 2/'^ 



f |2p6)|U)+ei=(i).2|y|2, 
m\'-^)+el={l)\V\^ 
m\^)+e^={l)\V\^ 

I |2p6)|tt)+eT-(|)-2|F|2, 



(15) 



- --\d^° 
2' 2/ 



|2/)| TT) + 

|2/)|i^; 



eT=2|T/|2, 
+ el^\V\\ 



(16) 



Multiplying by the transition rates to the intermediate 
state for all possible photon and electron spin polariza- 
tions. 



a+T = 




\\Ai?)\y\\ 


(17) 


= 


(2iAr 


v\\A,nv\\ 


(18) 


<y- T = 




^2\A,ny\\ 


(19) 


^- \ = 






(20) 



The total parallel and anti-parallel scattering is, 

13, 



TT= (fT+ T +0- i) 

ti= (a+ i +a- t) 
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\y\\ 
\y\\ 



(21) 
(22) 



leading to polarization (tt - ti)/(TT + Ti) = 
— (l/3)(5/6) for pure triplet states. 

The measured value of the polarization for each pho- 
toelectron energy gives an estimate of the amount of sin- 
glet and triplet character in the occupied states below 
Yjf- The experiment^'^^ find singlet character just be- 
low E p consistent with the t- J model and in contradiction 
to Ai holes that would Hund's rule triplet couple to the 
x^ — y^ electron. 

In our model, there are two types of Cu sites. The first 
is undoped with a single x"^ — y^ hole in a d^ state. The 
ejected photoelectron near the Fermi level comes from 
the Big combination of neighboring po- orbitals on the 
planar O that couples to the x^ — y^ electron in a singlet 
as described by Zhang and Ricci^ This is consistent with 
experiment and the t-J model. 

The second Cu is on a doped site with an out-of-plane 
polaron and a delocalized band comprised of x"^ — and 
Po- in our model. In this case, the final Cu d^ state has one 

and one x"^ — y^ hole with no spin correlation between 
them. Thus, tT=ti and the "polarization" arising from 
resonance scattering on doped Cu sites is zero. The only 
polarization observed arises from the undoped sites with 
singlet holes near the Fermi energy. 

The second experiment we consider is polarized xray 
absorption on La2-a;Sr2;Cu04 for x = 0.04 - 0.30,^2i^ A 
substantial O absorption with z-axis polarized xrays in- 
dicates there are holes in apical O Pz- In addition, xray 
absorption fine structure (XAFS)^!^ measurements ob- 
serve displacement of the apical O away from the Sr to- 
wards Cu consistent with hole formation in O Pz- Since 
the Pz hole character is compatible with our out-of-plane 
polaron assumption, we focus on the Cu result. 

The Cu absorption finds a few percent z^ character on 
the Cu sites. Our ab-initio calculations find the hole 



character to be approximately 85% of the x"^ 



y 



hole 
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TABLE 11: Polaron symmetries and energies from highest 
(most unstable electronic states) to lowest. The effective 
and Per orbital energies are taken to be 0. 



State 


Symmetry 


Energy 


P / P / 


E 






B2 


\J ^^ 4tp2 


S 


Ai 





P / P / 


E 


— \/2tp2 




B2 


J- /j-2 1 Ai2 
PP V PP ' ^^pz 


Not Polaron 


Aa 


Coupled to — 



character. It is too large compared with experiment. One 
could argue that the many-body response to the forma- 
tion of a Cu 2p core-hole is different for an undoped Cu 
versus a doped Cu where the delocalized x^—y^ band may 
suppress the white line due to the orthogonality catastro- 
phe or more strongly screen the core-hole potential. We 
are not convinced this is the sole reason for the small 
amount of hole character observed in the white line. 

A possible explanation is that the chiral polaron, 
spread out over four Cu sites as in figure ^ has more 
Pa and Pz character at the expense of from delocaliza- 
tion compared to the polaron centered around a single 
Cu site in figure [TJ A recent neutron pair distribution 
analysis^ is more compatible with a chiral plaquette po- 
laron. In this case, extracting a very small signal from a 
bulk average of bond distances and then using the mea- 
sured bond distances to infer orbital occupations is very 
model dependent. 



III. CHIRAL POLARONS 

The higher energy anti-bonding electronic states with 
apical Pz, , and p^ over a four Cu doped plaquette are 
shown in figure [31 The P^' and Pyi are degenerate. For 
simplicity, we have taken the two apical O Pz above and 
below each Cu and the Cu z"^ and 4s to be one Ai orbital. 
Thus, there are a total of eight states. The figure does 
not show the lower energy three bonding states (E and 
B2) since they are occupied. 

Table |TT] lists the energies of the eight polaron states 
for the case where the orbital energy of the "effective" 
Ai composed of and z^, is taken to be equal to the p^ 
orbital energy, Ez — ep — 0. There is an effective hopping 
matrix element, tpz, from pz to p^- tpp is the diagonal pa- 
matrix element. It is expected < tpz < tpp. 

Antiferromagnetic (AF) interaction of the polaron spin 
with the undoped Cu lattice renormalizes these cou- 
plings, but we expect P^;' and Pj,' to remain the most 
unstable electronic states. 

The effect of the undoped (fi spin background is seen in 
mean-field where the cP AF spins surrounding a plaquette 
are frozen with an "f spin on one sublattice and a | spin 
on the other sublattice. The additional energy of an S or 




FIG. 3: Projection onto CuOa of four polaron states and 
their symmetries. The higher energy anti-bonding states are 
shown. The fifth state of Aag symmetry is composed entirely 
of po- orbitals and is not a polaron state. This state becomes 
part of the delocalized — band at the unoccupied k state 
(tTjTt) and the occupied bonding — y^, pa band. We assume 
interaction with the undoped cP spin background, as shown 
in prior work,— i^i^i^i^ makes the chiral combinations P^i ±iPyi 
the most unstable electronic states (most stable hole states). 
When these chiral states percolate through the crystal, we 
further assume the electronic states composed of x'^ — y^ and 
Pct delocalize over these doped sites. There are a total of eight 
states. The lowest three have E {P^i and P^/) and B2 {Dxy) 
symmetry and are bonding combinations of the states in the 
figure. They are always occupied. The energies of the states 
are listed in table lU 

Dxy polaron with spin a due to AF coupling of a with 
the cP spins is zero since the average <P spin seen by the 
polaron is zero. For P states, the polaron spin couples to 
one sublattice. a can be aligned with the sublattice spin 
leading to a further destabilization of the P state. 

P hole states were found in the exact diagonalization 
of Gooding^ for a t-J model on a 4 x 4 lattice with 
an additional hole allowed to delocalize on the interior 
2x2 lattice. This is in accordance with theoretical 
predictions 

Based on the energies in table [TTl the mean-field de- 
scription of the (P spins, and exact results on a 4 x 4 
lattice^ we assume the polaron hole has P symmetry. 

For a single hole delocalized in a small region of an AF 
spin background, it has been showit^i^ii the chiral states 
\/^/2{Px' ± iPyi) are the correct spontaneous symmetry 
breaking states for the hole rather than P^^' and Pj,/ , be- 
cause the complex linear combinations are compatible 
with the long range twisting of the AF lattice spins into 
a stable configuration topologically distinct from the AF 
ground statei^ 

In this paper, we assume doping introduces hole char- 
acter in out-of-plane orbitals that can delocalize over a 
small number of sites in the vicinity of the dopant. The 
most favorable configuration for the polaron is taken to 



8 



be the chiral state. If there was a single dopant in an in- 
finite (P crystal, then the chiral states, l/^/2{Px' ±iPyi), 
would be degenerate. These two states are time-reversed 
partners. 

In a finitely doped system, the environment of each 
polaron is different and the two chiral states may have 
different energies. We assume, in a doped cuprate, the 
chiral states are the correct polaron eigenstates, but the 
energies of the two states may be different. This leads to 
a model of the polarons where the splitting between the 
chiral states along with all the other states represented 
in table [Tll and figure [3] are distributed differently for 
each plaquette. The assumption of a completely uniform 
probability distribution of different polaron state energies 
throughout the crystal leads to neutron uj/T scaling as 
shown in section fVAl A linear resistivity, derived in sec- 
tion IVIIi arising from the Coulomb scattering of — y"^ 
band electrons with the polarons is also obtained with a 
uniform energy distribution. 

This model of non-degenerate chiral polarons implies 
time-reversal symmetry is broken. At any instant, the 
number of "up" chiral polarons should equal the num- 
ber of "down" chiral polarons and macroscopically the 
cuprate is time- reversal invariant. There is recent exper- 
imental evidence for local time-reversal symmetry break- 
ing in neutron scatteringi^ 



IV. PERCOLATION 

There are three basic assumptions of our model. First, 
doping leads to additional holes in out-of-plane orbitals 
that form chiral states as shown in figures [2] and [S] 

Second, when these polaron plaquettes percolate 
through the crystal, a band is formed with the x'^ — 
and po- orbitals on the percolating swath. This metallic 
band interacts with the — hole spins on the un- 
doped Cu sites and the plaquette polarons. The random 
distribution of impurities leads to a distribution of the 
energy separations of polaron states shown in figure [H 

Third, this energy distribution is uniform. The linear 
resistivity arises from this assumption as shown in section 
IVIII Since the resistivity is non-linear for certain dopings 
and temperature ranges, this assumption is not always 
valid. 

The transition from spin-glass to superconductor in 
La2-2;Sr2;Cu04 at a; sa 0.05^" and from an antiferromag- 
net to superconductor at (5 « 0.35^ in YBa2Cu306-i-5 
occur at the doping when the polarons percolate through 
the crystal. 

In this section, the site percolation doping values are 
computed for La2-xSra;Cu04 and YBa2Cu306+5. Rea- 
sonable assumptions for the distribution of plaquettes 
are used to approximately simulate the repulsion of the 
dopants. The computed values are close to known phase 
transitions in these materials. 

We also computed the percolation for two additional 
systems. The first is La2-2;Sra;Cu04 where each Sr dopes 



exactly one Cu site as shown in figure [T] and the second 
is a 2D square lattice with plaquette doping. The com- 
puted La2-a;Sr2;Cu04 1-Cu percolation value of x w 0.20 
is associated with the observed orthorhombic to tetrago- 
nal phase transition.^- 

For the 2D square lattice with four Cu plaquette dop- 
ing, percolation occurs at x w 0.15. We believe the 
2D percolation of the plaquettes should be associated 
with the transition from insulator to metal at x « 0.15 
found by low temperature resistivity measurements in 
large pulsed magnetic fields. This is further discussed 
in section IVlTl 

All percolation computations described here were per- 
formed using the linear scaling algorithm of Newman and 
Ziflf.^o 

In all these calculations, we simplify the problem 
by using Cu sites only. For La2-a;Sra;Cu04, we take 
each Cu to have four neighbors in the plane at vectors 
(±a, 0, 0), (0, ±6, 0) and eight neighbors out of the plane 
at (±a/2, ±&/2, ±c/2) where a, fe, and c are the cell di- 
mensions. Thus, each Cu has a total of 12 neighbors in 
the site percolation calculations. 

For all YBa2Cu306+5 calculations, we take each pla- 
nar Cu to be connected to a total of six Cu atoms. There 
are four nearest neighbors in the same Cu02 plane, one 
neighboring Cu on the adjacent Cu02 across the inter- 
vening Y atom and one Cu on the neighboring chain. 
The chain Cu is connected to the two Cu atoms in the 
Cu02 planes above and below itself. We assume a chain 
O dopes three Cu atoms, two in Cu02 planes and the 
corresponding chain Cu as shown in the constraints of 
figure m 

Table IIIII lists the computed percolation values for 
La2-2;Sr2;Cu04, a 2D square lattice, and YBa2Cu306+5 
for various types of doping and doping constraints. These 
constraints were chosen to simulate the repulsion of the 
dopants and are approximations to the actual distribu- 
tion of dopants in the cuprates. 

The La2-KSr2;Cu04 and YBa2Cu306-(-5 calculations 
are on 200 x 200 x 200 lattices with 2, 500 different dop- 
ings. The square lattice size is 2, 000 x 2, 000 with 5, 000 
different dopings. 

The first La2-a;Sra;Cu04 calculation is the critical dop- 
ing for percolation of doped Cu where each Sr dopes the 
single Cu shown in figure [1] instead of the four Cu pla- 
quette of figure [2l Although we assume plaquettes are 
created at low dopings, once the doping is large enough, 
there is crowding of the plaquettes. Single Cu polarons 
are formed. This single Cu percolation calculation is 
an approximate measure of the doping for the transi- 
tion from predominantly doped plaquettes to single site 
polarons. A phase transition at this crossover doping is 
expected. The computed percolation of x w 0.20 matches 
the orthorhombic to tetragonal transitioniS doping. 

From the table, the critical doping for 3D plaquette 
percolation in La2-a;Sra;Cu04 is x « 0.05 regardless of 
the applied doping constraints and matches the spin-glass 
to superconductor transitions^ 



TABLE III: Percolation values for various structures and dop- 
ing scenarios. The constraints column references the fig- 
ures describing the applied constraint. All La2-a;Sra;Cu04 
and YBa2Cu306-i-i results are obtained for a 200 x 200 x 
200 lattice with 2, 500 doped ensembles. The percolation 
value is the computed critical x in La2-a;Sr2,Cu04 and S in 
YBa2Cu306+i. The square lattice results are for a 2,000 x 
2, 000 lattice with 5, 000 ensembles. The digit in parentheses 
is the error in the last digit. 



Structure 


Dopant Type 


Constraints 


Percolation 


LSCO 


1-Cu 


none 


0.19617(2) 


LSCO 


4-Cu 


none 


0.05164(1) 


LSCO 


4-Cu 


fig. Hi 


0.05097(1) 


LSCO 


4-Cu 


fig. |li,b 


0.04834(1) 


LSCO 


4-Cu 


fig. |4li,b,c 


0.04880(1) 


LSCO 


4-Cu 


fig. |4^,b,c,e 


0.04904(1) 


LSCO 


4-Cu 


fig. |lJi,b,e 


0.04862(1) 


LSCO 


4-Cu 


fig. |4li,b,c,d 


0.04926(1) 


LSCO 


4-Cu 


fig. |4^,b,c,d,e 


0.04943(1) 


Square 


4-Cu 


fig. |li,b 


0.15053(1) 


YBCO 


3-Cu 


fig.[5K 


0.31162(2) 


YBCO 


3-Cu 


fig-Et 


0.32890(2) 


YBCO 


3-Cu 


fig. [S1d,c 


0.36098(2) 



This is because the plaquette percolation values are 
approximately 1/4 of the single Cu percolation result of 
X « 0.20. 

For YBa2Cu306+5, the more realistic doping con- 
straints are the second and third cases where 6 « 0.33 
and S ~ 0.36 since O chains should not have a preference 
for which Cu triple to dope. Experimenl^ finds S « 0.35. 

From these results, we conclude the plaquette po- 
laron model with percolation can obtain known insu- 
lator to metal phase transitions in La2-2;Sr2;Cu04 and 
YBaaCugOe+i. 



V. NEUTRON SCALING AND 
INCOMMENSURABILITY 

A. Scaling 

Neutron spin scattering measures the imaginary part 
of the magnetic susceptibility, x(q, i^)- 

The integral of the imaginary part of the spin sus- 
ceptibility /d^q x"(q, 1^) over the Brillouin zone where 
X = x' + ix" has been foundi8a9-20,2is^^ to be a 
function of u/T. The integral is the on-site magnetic 
spin susceptibility. 

The oj/T scaling is unusual because x" ^ lu/ Jdd for an 
antiferromagnet and x" ^ Lo/Ep for a band where Jdd 
is the AF spin coupling and Ep is the — band 
Fermi energy. 

In this section, we show that the single polaron suscep- 
tibility is a function of lo/T when the energy difference 
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FIG. 4: Applied Sr doping constraints used for the 
La2-iSr3:Cu04 and 2D square lattice plaquette percolation 
calculations shown in table IIIII For each four Cu plaquette 
in figures [2] and |3l two Sr atoms, one above and one below 
the plaquette (upper and lower Sr), can dope the Cu's. Each 
figure shows the disallowed configuration of Sr doping. It is 
assumed all ±90° and 180° rotated configurations are equiva- 
lent to the figure and also disallowed, (a) an upper and lower 
Sr doping the same plaquette. (b) a Cu atom in a plaque- 
tte doped by two different Sr atoms. This figure includes the 
three cases of two upper Sr, two lower Sr, and one upper and 
one lower Sr. This is the no overlap constraint, (c) adjacent 
plaquettes doped by two upper Sr or two lower Sr. (d) ad- 
jacent plaquettes doped by one upper Sr and one lower Sr. 
(e) nearest neighbor upper Sr and lower Sr in different LaO 
planes. 



between polaron chiral states with opposite spins and 
chiralities is uniformly distributed. 

The q dependence of the polaron susceptibility, 
Xp (q, w), is peaked at q = if spin-flip polaron scattering 
dominates at low energy. x!p{^.i ^) is peaked at q = (tt, tt) 
if the polaron spin and chirality flip at low energy. The 
latter scatters the polaron to its time-reversed state. The 
time-reversed chiral polarons are the low energy excita- 
tions, as shown in section IV Bl The q dependence is 
broad because the polaron is localized over a four-site 
plaquette. Since the total susceptibility is dominated by 
q near (tt, tt), the polaron susceptibility is approximately 
momentum independent, Xp(ctii^) ~ 'x'pi'^)- 

The imaginary part of the polaron susceptibility is 
found to be of the form x'pi^) ^ tanh(/3a;/2) and sat- 
isfies scaling. The coupling of the polaron spin and chi- 
ral orbital state to the undoped (P spins causes the total 
susceptibility to become incommensurate. This is shown 
in the next subsection IV Bl where we compute the static 
spin structure factor for classical spins perturbed by chi- 
ral polarons. 

In this section, we show that coupling to the undoped 
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FIG. 5: Doping constraints used for YBa2Cu306+5 calcula- 
tions shown in table IIIII The chain O and its three doped 
Cu sites are shaded gray. Three diflterent doping scenarios 
are shown, (a) the chain O dopes a fixed Cu triple, (b) the 
chain O randomly dopes one of the two possible Cu triples. 
No triple may be doped by two adjacent chain O. (c) If two 
chain O are in the same cell, the two Cu doping configurations 
where the triples are closest to each other are not permitted. 



spins leads to a dynamic susceptibility consistent with 
the measured form^*''^^ in equation [551 

Consider a polaron as in figure [S] with energy A sepa- 
rating the down spin state from the up spin and occupa- 
tions rij and ?i| in thermal equilibrium. Since the polaron 



is always singly occupied, + 
Solving for and n^, 



1 and n 



p/3A 



71^ = /(-A), nT=/(A), 
where /(e) is the Fcrmi-Dirac function, 

m ' 



(23) 



(24) 



To compute the dynamical polaron spin susceptibil- 
ity, consider an applied magnetic field _Be*"* normal to 
the spin quantization axis. The alternating field induces 
transitions between the two states if A = uj. 



m 

ff-AJ 



FIG. 6: Energy difference A between an up and down spin 
polaron. /(—A) and /(A) are the equilibrium occupations of 
the spins states where /(e) is the Fermi-Dirac function, /(e) = 
l/(e'''^ -I- 1). Beff is the effective magnetic field that splits 
the spin energies by A. The chiralities of the two states are 
different for low energies leading to a broad peak in Xp(q,'^) 
at (tt, tt). 



Let W — = be the induced transition rate 
between the two spin states. The total absorption rate 
is, 

{P{^))^h^^W{n^-n^). (25) 
Using equation [23l 

(P(A))=a^M^[/(-|A|)-/(|A|)], (26) 



(P(A)) = ;iwVKtanh(/3|A|/2). 



(27) 



The absolute values of A are used above because the 
absorption rate is independent of which spin state is lower 
in energy. The transition rate W is. 



(28) 



Averaging W over all spin quantization directions multi- 
plies equation [28] by 2/3, 



(29) 



Let p(A) to be the probability distribution of energy dif- 
ferences for spin and chirality flips. Summing over all 
polarons, the total absorption rate is, 

(P(w)) =htj [ dA p(A)tanh(/3|A|/2)(W^(A)), (30) 



{P{cj)) = [ ^) u;^,lB' [p{u;) + p{-u;)] tanh 



(P(c.)) = ('|)^^^i?V(c^)tanh(^ 



(31) 
(32) 



where p(— A) = p(A) because there is an equal number 
of polarons with an up spin lower in energy than a down 
spin as there are down spins lower in energy than up 
spins. 
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The absorption rate can be written in terms of the 
imaginary part of the polaron susceptibility as, 



(33) 



leading to the imaginary susceptibility per polaron of, 



x;'M-(^)A^|pHtanh(^). (34) 



The probability density of polaron energy separations 
p{A) is taken to be uniform with J c?Ap(A) = 1 and of 
the form, 



p(A) = 



— J 2A„ 



0, 



|A| > A,- 



(35) 



where A^ax is doping dependent, A^ax = Amax(a;). 

Equations [34] and [35] show that the polaron suscepti- 
bility is a function olu/T and has the approximate form 
seen in experiments.— The functional form of x^' in- 
creasing from Xp (0) = and saturating for j3uj > 1 arises 
from the thermal occupations of polaron states with en- 
ergy splitting u!. When Puj <C 1, the two states have al- 
most equal occupation and the absorbed energy is small 
from equations and For /3uj ^ 1, the lower en- 
ergy state is always occupied and the higher energy spin 
state is always unoccupied. In this case, the absorption 
saturates. Since the polaron spin density p(A) is con- 
stant up to Amax, it is Plu that determines the amount of 
absorption due to the difference of the two Fermi-Dirac 
occupation factors. Finally, if there are no spin flip ener- 
gies smaller than Amin, then Xp{^) is zero for uj < Amin- 

The measured susceptibility for La2-a;Sr2;Cu04 at a; = 
0.04 is normalized and fit by the expressio n^°i^^ 



tan-i[ai(/3w) +a3(/3a;)3]. 



(36) 



with ai = 0.43 and 03 = 10.5. This curve rises to the 
saturating value of one faster than our expression in equa- 
tion [33 

The contribution to the susceptibility from the un- 
doped dP spins and the metallic — y'^ band has not 
been included. The band contribution is on the order 
of Lu/Ef where Ef is the Fermi energy and can be ne- 
glected. The imaginary susceptibility from the undoped 
cP spins is on the order™ of cu/T where F is several Jdd 
to Ef and can also be neglected. The real part of the (P 
susceptibility is approximately constant up to the energy 
w ~ F. We may therefore take the (P susceptibility to be 
real and uj independent for small uj. Thus, the uj depen- 
dence of the total susceptibility arises from the polaron 
susceptibility in equation 1341 

The q dependence of the susceptibility is incommensu- 
rate from the calculations of the next section and of the 
formjii 



x(q) = 



i + (q- 



Qf(e/«)" 



(37) 



+ 




0^ 



+ 



u 



FIG. 7: Diagrams to sum in the random phase approximation 
for the dynamic susceptibility x(ci)<^)- The first term is the 
static susceptibility x(ci) in equation [33 and the shaded loop 
is the polaron susceptibility X^(io) with imaginary part shown 
in equation 1341 U is the coupling energy. 



where Q is the incommensurability peak vector and ^ 
is the correlation length. From the computations in 
the next section, Q is shifted from (7r/a,7r/a) along 
the Cu— O bond direction and agrees with experiment. 
^ « a/ y/x is the mean separation between Sr. 

Summing the random phase approximation diagrams 
in figure [7] leads to. 



"x(q,^^)' 


-1 


[X(q)] 


-1 




. /4 . 




. /4 . 







(38) 



x(q,^) 



l + (q-Q)'(e/a)' 



XoU'^ [Xp(w)/Ai 



where we have used (^/a)^a; = 1 for the Xp term. 
Using the integral. 



(39) 



Im 



A + Bq^- iC 



tan 



(27rC) qdq 
{A + Bq-^f + C^ 
C 
A 



( — 1 tan 

\B 



(40) 



and defining Xp(uj) = Xp(w)//i^, the integrated imagi- 
nary susceptibility is. 



J X"{<i,i^)dq = 7r/x|xo tan" 



1 - XoU^xU^} 



(41) 



In the above expression, Xp = Xp+*Xp been expanded 
into real and imaginary parts. 

Using x'pi-^) = X'pi^) and Xpi-^) = -Xp{^), the 
equation can fit the experimental curve in equation 1361 



B. Incommensurability 

There is an antiferromagnetic coupling, Jpd, between 
the polaron spin Sp and the neighboring d^ spins. The 
polaron is dclocalized over four Cu sites. The probabil- 
ity of the hole residing on a particular Cu is 1/4 lead- 
ing to the estimate Jpd ~ (1/4) Jdd where Jdd is the 
undoped d^ AF coupling. The effective coupling of a 
chiral polaron to the d^ spin background is knowniiSiii^ 
to induce a twist in the neighboring spins. This can be 
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FIG. 8: Coupling of chiral polarons to neighboring spins. 
Sp is the polaron spin and it couples to four pairs in the 
cyclic order shown by the arrows and equation 1421 The chiral 
coupling for the P^i — iPyi polaron reverses the cyclic order- 
ing of the spins as seen in the figure and leads to the same 
expression as the P^i + iPyi with J^h — >■ —Jch- 



encapsulated in a topological charge term^ of the form 
=FJch[Sp • (Si X S2)] where Sp is the polaron spin and 
Sdi, Sd2 are adjacent (P spins as shown in figure [51 

The expectation value (Sp • (Si x S2)) = for states 
invariant under time reversal. Thus, the expectation 
value of the topological charge is zero for the real polaron 
states Pa;', Pyi , S, and Dxy The complex linear combi- 
nations in the chiral states lead to non-zero topological 
charge. The above chiral coupling term is the simplest 
coupling of chiral polarons to the neighboring spins. 

The referencesi^ii^ considered holes in both the t-J 
and three-band Hubbard models that can delocalize over 
a four-site plaquette. In our model, — y^ spins delocal- 
ize on the plaquettes forming a band when the polarons 
percolate. Our chiral coupling is between a polaron spin 
and the adjacent spin sites that may be undoped dP or 
another polaron spin. The specific form for the coupling 
is analogous to previous work. 

The coupling of the x'^ — band to the neighbor- 
ing spins is smaller than the chiral coupling of the po- 
laron and (P spins. The perturbation arising from the 
band spin coupling selects incommensurability along the 
Cu— O bond directions as shown at the end of this sec- 
tion. 

The polaron chiral coupling of P^' ± iPy' to the neigh- 
boring cP spins is, 

Hch = TJch {Sp • (Si A S2) + Sp • (S3 A S4)+ ,^2) 
Sp-(S5AS6) + Sp.(SrAS8)}, 

where Jch > and the spins are labelled in figure [51 The 
antiferromagnetic coupling of polaron spins to undoped 
(f spins is 

Hpd - JpdSp • (Si + . . . + Sg) . (43) 

Electronic hopping matrix elements are on the order of 
0.5 — 1.0 eV. The chiral coupling, Jch, is estimated to be 
less than or of the same order. Gooding et alfS obtain 



Jch ~ SJdd from numerical simulations and by comput- 
ing the effective next-nearest neighbor antiferromagnetic 
coupling, J^^, induced by chiral polarons at very low dop- 
ing. J^^ is then compared with Raman data to obtain 

Jch- 

In the spin Hamiltonian of Gooding et al.^^ the chi- 
ral coupling term is squared, — Jch[Si • (S2 A Sa)]^, in 
contrast to our linear terms in equation 1421 The over- 
all energy scale of Jch is similar. We take Jch = 3Jdd 
where Jdd ~ 0.1 eV in our computations of the static 
neutron structure factor. We have found our results for 
the magnitude of the incommensurability are indepen- 
dent of the precise values of all of the parameters. The 
only necessary feature to obtain incommensurability is 
that the chiral coupling, Jch, is sufficiently large to break 
the (tt, tt) spin ordering from the antiferromagnetic spin 
coupfing, Jdd- 

There is an antiferromagnetic spin-spin coupling be- 
tween the polarons. Jpp is the coupling between Spi 
and Sp2 and Jpp is the coupling between Spi and Sp3 
shown in figure [SI An estimate of Jpp and J^^ is ob- 
tained in a similar manner to Jpd- For Jpp, the polarons 
have two adjacent pair sites. An antiferromagnetic cou- 
pling occurs for every adjacent pair. This occurs with 
probability (1/4)^ = 1/16. There are two pairs for Jpp 
and one for J^^ leading to estimates Jpp « (1/8) J^d and 
J^p « (1/16) Jdd. 

Figure [9] shows various chiral couplings when polarons 
are adjacent to each other. Using a similar analysis, we 
estimate the magnitude of the chiral couplings to be, 
Jppp ~ (l/4)Vc/i, Jppd ~ (1/4) Jch, and J^^^ « (1/4) J^^. 

The total spin Hamiltonian for the spins and po- 
larons is 

H^Hdd + Hpd + Hpp + Hll' , (44) 

where Hdd is the antiferromagnetic spin-spin coupling 
with Jdd ~ 0.1 eV. Hpd is the polaron-d^ coupling and 
Hpp is the polaron-polaron spin coupling. is the 

total chiral coupling. 

The chiral coupling iJ*^* is invariant under polaron 
time reversal that fiips the chirality of a single polaron 
Px' ± iPy' Px' T iPy' or {Jch — > - Jch) and the po- 
laron spin, Sp — Sp. Hpd — > —Hpd is not invari- 
ant under time reversal of the polaron. When the chi- 
ral coupling is much larger than all spin-spin couplings, 
Jch > Jdd > Jpd > Jpp > Jpp, the ground state energy 
becomes independent of the polaron chiralities. 

This is an important point because it means the energy 
to simultaneously flip the chirality and spin of a polaron 
has an energy scale Jdd while flipping either the chirality 
or the spin, but not both, has an energy scale Jch- 

Figure fTUI shows the minimized spin ordering surround- 
ing a polaron when the chiral coupling, Jch, dominates 
the polaron spin to coupling, Jpd- Increasing Jch fur- 
ther does not change the spin ordering. This leads to an 
incommensurability that is weakly parameter dependent. 
All the neighboring d^ spins in the figure are orthogonal 
to the polaron spin. The Hpd antiferromagnetic energy 
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FIG. 9: Schematic of adjacent polaron configurations. Jppp is 
the chiral coupHng between three polarons, Jppd couples two 
polarons to a spin, and J'pp^ couples two polarons to a 
spin when one polaron is shifted by one lattice spacing. The 
couplings are shown for the case where all the polarons are 
P^i + iPyi . For each opposite chirality polaron, P^i — iPyi , 
the coupling should be multiplied by —1. The figure does not 
exhaust all possible couplings. For example, there is another 
term Jppp if Sp4 is shifted to the right by one lattice spacing. 




FIG. 10: Minimized energy configuration for a polaron with 
spin pointing out of the page surrounded by 8 cfi spins. Only 
Jdd and Jch are non-zero. This represents the regime domi- 
nated by Jch- 

is zero. The energy difference of the time reversed po- 
laron in the same background is also zero. The spin-spin 
couplings, Jpd etc, lead to non-zero energy differences. 

The static neutron spin structure factor is computed 
by minimizing the energy in equation [44] on a finite 2D 
lattice with classical spins of unit length, — 1. 
Each undoped (P site has a spin and every polaron has 
an orbital chirality, Px' ± iPy', and spin. All the terms 



for the classical Hamiltonian are described above along 
with the parameters used. The only constraint on the 
polarons is that they may not overlap, but otherwise they 
are randomly placed in the lattice. At this point, the 
effect of the delocalized — y^ band electrons is ignored. 

One can imagine additional constraints on the place- 
ment of the polarons arising from polaron-polaron 
Coulomb repulsion. Also, calculations on 3D lattices 
with a small interlayer antiferromagnetic spin coupling 
can be done. The addition of polaron constraints and 
the third dimension does not to change the computed 
incommensurability or its location in the Brillouin zone. 
These effects are ignored in this paper. 

Finally, calculations with periodic and non-periodic 
boundary conditions were performed to ensure there is 
no long range twist in the spins that is frustrated by 
periodic boundary conditions. No major difference was 
found in the computed structure factors. This is likely 
due to the small energy difference between a polaron and 
its time reversed partner. 

Computations were done on 256 x 256 lattices with 
polaron dopings of x = 0.075, 0.10, and 0.125. A ran- 
dom configuration of polarons was chosen subject to the 
constraint that no polarons overlap. Starting spins and 
chiralities are randomly generated and the initial energy 
is calculated. The energy is minimized by performing 
local minimizationSf^^ 

A spin is selected and the effective magnetic field on the 
site is computed. Since the Hamiltonian in equation 031 
is linear in the spins, the energy arising from the chosen 
spin is minimized by aligning it with the magnetic field. 
If the spin is a polaron spin, then the effective magnetic 
field is computed for both orbital chiralities to determine 
the chirality and spin that minimizes the energy. The 
chirality is flipped if a lower energy can be obtained. The 
program loops through all the spins, determines the new 
energy, and compares it to the previous energy to decide 
on convergence. Calculations were performed on 5,000 
different polaron configurations for each doping value. 

The Grempel algorithn>i^iii was used to obtain the 
global minimum. This algorithm is similar to raising the 
temperature to allow the energy to climb out of local 
minima and then annealing. Unlike Gooding et al.j^ we 
found the Grempel steps lowered the energy minimally 
and made no difference to the static neutron structure 
factor. It is likely that this is due to our linear chiral 
coupling in equation 1421 A squared chiral coupling, used 
by Gooding et al., makes the minimization more diffi- 
cult and computationally expensive. Thus, we were able 
to minimize larger lattices and more ensembles to obtain 
smaller error bars on the results. Our calculations con- 
stitute a different physical model than Gooding et alJ^ 
despite the similar computational methods used. 

Finally, we found that including the polaron-polaron 
spin and chiral couplings shown in figure [9] does not alter 
the results. The dominant couplings in terms of the min- 
imized spin structure are the (P spin coupling Jdd and 
the chiral coupling Jch- The results shown below exclude 
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FIG. 11: Static spin structure factor for spins at doping 
X — 0.10. The structure factor is incommensurate with a ring 
of radius k m {211 /a)x around {n / a,n / a). The total sum over 
the Brillouin zone satisfies the normalization, A''"^ S{k) — 
1 where A'' is the total number of cells, (a) 3D plot of structure 
factor, (b) contour plot centered at (ji /a,-n /a). 



any chiral couplings involving more than one polaron. 

Figures [TT] and [T2l show our results for the spin correla- 
tion at dopings x = 0.075, 0.10, and 0.125 on a 256 x 256 
lattice averaged over an ensemble of 5, 000 configurations 
for each doping. The structure factor is dimensionless 
and is normalized such that its integral over the Bril- 
louin zone is one, S{k) = 1, where N is the to- 
tal number of cells. Figure [T51 shows part of a minimized 
spin structure at x = 0.10. 

Due to the large number of ensembles, the error bars 
for the plotted values are small. For x = 0.075, they are 
±0.12, ±0.15, and ±0.12 for the diagonal peak, (7r,7r), 
and the Cu— O bond peak respectively. For x = 0.10, 
the error is ±0.09, ±0.08, and ±0.07. For x = 0.125, the 



error is ±0.07, ±0.05, and ±0.05. The error decreases 
sharply and becomes negligible on the scale of the figure 
as k moves past the peaks and farther from (7r,7r). 

From figure 1121 the diagonal incommensurate peak 
is shifted from (7r/a,7r/a) by (27r/a)(a;/\/2)(l, 1) and is 
of length {2TT/a)x. The peak along the Cu— O bond 
direction is shifted in the range (27r/a)a;(l/V2, 0) to 
(27r/a)x(l, 0). The Cu— O bond direction shift is exper- 
imentally observed for the metallic range x > 0.05 ^°i^'^ 
in La2_a;Sra;Cu04 and the diagonal shift is seen in the 
spin-glass regime 0.02 < a; < 0.05.^'*: 

Since the difference between (27r/a)a:(l/-\/2, 0) and 
(27r/a)a;(l, 0) is small, it is difficult to resolve the precise 
peak doping value within our finite size computations. If 
the structure factor derives from broadened Lorenztians 
centered at the four diagonal points around (7r/a,7r/a), 
then a peak in the Cu— O bond direction would be ex- 
pected at (27r/a)a:(l/\/2, 0) from the two closest peaks. 
The contributions from the remaining two peaks shift the 
peak closer to (7r/a,7r/a) rather than away from it as is 
seen in figure 1121 Thus, the best we can currently say 
with the calculations is there is a ring of incommensu- 
rate peaks approximately a distance k ~ (27r/a)x from 
(tt/o, ir/a). 

From the widths of the peaks, the correlation length 
is approximately the mean separation between the po- 
lar ons, a/y/x. 

We present a heuristic derivation for why the spin 
structure factor is incommensurate with a shift from 
(7r/a,7r/a) of magnitude (2TT/a)x. Similar to previous 
workfi^ our calculations find that the minimum spin con- 
figuration consists of undoped patches of (P spins aligned 
antiferromagnetically with the polarons acting to rotate 
the direction of the antiferromagnetic alignment of adja- 
cent patches. This is seen in figure [TS] 

Consider an area A. The number of polarons in this 
area is Np = Ax/a^. When the chiral coupling domi- 
nates, the effect of a single polaron on the neighboring 
(P spins is shown in figure [TUl The polaron rotates each 
adjacent spin on opposite sides of the polaron by tt, or 
27r total. If this net 27r rotation is rigidly transmitted 
to an antiferromagnetic patch, then the polaron rotates 
a patch by an angle, 59 = 27r/(l/x) = 27ra;. Thus, the 
net rotation per spin is NpS6. If the area A is chosen 
such that the net rotation per spin is 27r, or NpS9 — 27r, 
then A — {a/x)"^. A translation by L = ^/A = a/x re- 
turns to an identical antiferromagnetic patch. This leads 
to a shift of the spin correlation peak from (7r/a,7r/a) to 
5Q = 27r/L or 5Q = (27r/a)a;. 

C. Kinetic Energy of the — y'^ Band 

The energy contribution from the delocalized x^ — 
band electrons has not been included in the minimiza- 
tion of equation 1441 A complete minimization would also 
compute the change in the band energy to determine the 
direction to align a given spin during our sweep through 
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FIG. 12: Static spin structure factor for d spins at dop- 
ings X = 0.075, 0.10, and 0.125 in La2-^Sra;Cu04. Each 
figure plots S'(fc) starting from k = (0,0) to {Ti/a,-n/a) and 
then to (0,7r/a). The structure factor is incommensurate 
with a ring of radius k ~ (2-K/a)x around (ji / a,n / a). The 
vertical lines are drawn to highlight specific incommensu- 
rate vectors. If the structure factor was purely derived from 
the sum of four Lorenztians along the diagonals a distance 
{2Ti/a)x from {n/a,n/a), then the peak along the Cu— O 
bond direction would be slightly less than the shift shown at 
SQ = (27r/\/2t[, 0). The normalization is the same as figure 



the lattice spins. This effect is included in mean- field 
below. 

A dP spin ordering of momentum q hybridizes band 
electrons of momentum k and k -I- q with a coupling en- 
ergy V on the order of Jdd ~ 0.1 eV. This mixing of k and 
k -|- q perturbs the band energies and the total ground 



state energy. We calculate the band energy change in 
La2-2;Sr2;Cu04 for a ring of q vectors at the computed 
incommensurate length q = {2iT/a)x. The q vector pro- 
ducing the lowest band energy is the observed neutron 
incommensurability. 

The idea that the band kinetic energy change in the <P 
spin background determines the final neutron incommen- 
surability has been suggested by Sushkov et al.^'"^'^^'^^ for 
the t — t' — t" — J model. In their model, the magnitude 
of the incommensurability arises from the doping x, the 
antiferromagnetic cP spin stiffness Ps, the hopping ma- 
trix element t, and the quasiparticle renormalization Z 
where q = {Zt/ ps){'2,'K / a)x. Their self-consistent Born 
approximation calculation of the quasiparticle dispersion 
in a spin-wave theory background finds values for the 
parameters such that Ztj ps ~ 1. The magnitude of the 
incommensurability is less dependent on the detailed pa- 
rameters for our Hamiltonian in equation 1441 

For a given spin incommensurability vector q and to 
— if' band coupling, ~ 0.1 eV, the Green's function 
satisfies, 



G-i(k,t^) = Goi(k,c^)-y2^o(k + q,c^) 

- l/2Go(k-q,u;). (45) 

The vector — q must be included with q because the cou- 
pling Hamiltonian is Hermitean. Solving for G(k,a;), 

G(k,c^) = Go(k,c^){l-^2go(k,c^) 

X [Go(k + q, c^) + Go(k - q, uj)]}-^ . (46) 

Expanding to order V'^, 

G{k,u;) = Go(k,c^) + F2Go(k,c^){Go(k + q,c^) 
+ Go(k-q,w)}Go(k,w) 
+ 0{V^). (47) 

The number of electrons in the k state up to energy e is, 



n(k, e) 



duj 



-liR G(k,tj) 

TT 



(48) 



X lm{Go(k,cj) [Go(k + q,w) 

+ Go(k-q,c^)]Go(k,c^)}, (49) 

where no(k, e) is the unperturbed occupation. The total 
density of states per spin at energy e is. 



k ^ ^ k 

A percolating band has Green's function, 

1 



Go(k,w) 



U! - ek + iTk 



, rfc>o, (51) 
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FIG. 13: Projection of spins onto xy plane for an a; = 0.10 minimized structure. The spin staggered magnetization is along 
the z-axis out of the plane of the paper and the magnetization of the spins is along the x-axis. Only a 30 x 30 subset of the 
256 X 256 lattice is shown. The undoped spins are shown in red and the polaron spins are green. The four Cu sites of each 
polaron are indicated by black dots. 



where is the linewidth. Using equations H51 and 

«o(k,c.) = i + ltan-i(^^^^). (52) 

The perturbation shifts the Fermi level to e/ + 5e/. 
The total number of electrons is conserved, leading to 
the equation for 6ef, 

k k 

The total energy of the band electrons per spin is given 
by, 

^tot(q,'^^)- /'^ 'dc^c^^«(k,c^). (54) 

J-co ^ 

In appendix [Xl the expressions for Sef and Etoti^i V) to 
order are derived, 

6efNo{ef) + V^ H ^ [i?(k, k + q, c^) 



i?(k,k-q,w)] =0, (55) 



Etot{c\, V)=Eg + V^ dLo{u - ef) 



i?(k, k + q, w) + R{k, k - q, cj) 



(56) 



where No{ef) is the unperturbed density of states per 
spin, Eg is the unperturbed {V = 0) energy, and 



i?(k,p,(^) 



Im 



{uj - ek + iTkY 



(w - + iTp) 



(57) 



The integral /I^ -R(k, p, u;)duj can be evaluated an- 
alytically, thereby allowing us to accurately compute 
the small energy change of the Fermi energy 5ef and 
Etotil, V). This is done in appendix \K\ 
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The band energy is given by, 



gfc = — 2ii(cos fc^ + cos /c, 
—2^2(003 2k 



,y) — Atli COS cos ky 
cos2ky) + ex2_y2. 



(58) 



We use the band structure parameter o 1 ' for 
La2-^Sr^Cu04 given by ti = 0.25 eV, tu = -0.025 eV, 
and t2 = 0.025 eV. ti is the nearest neighbor hopping, 
til is the next-nearest neighbor diagonal term, and ^2 is 
the hopping along the Cu— O bond direction from two 



0.10, 



0.133 eV leads 



lattices site away. For x 
to e/ = 0. 

We calculated the Fermi energy shift and energy for 
incommensurability along the diagonal and Cu— O bond 
direction of magnitude q = {2'K/a)x a,t x — 0.10 with 
V = 0.1 eV. The electron linewidth, Ffc, is chosen to be 
the sum of an s-wave and d-wave term, 



Fs + T£i{coskx — cos kyY, 



(59) 



where Fs = 0.01 eV and F^ = 0.01 eV. 

The addition of a d-wave term to the linewidth arises 
from the k dependence of the Coulomb scattering rate 
with polarons discussed in section IVII Al 

The energy changes and Fermi level shifts in eV are, 



= -0.01456, (5e/ = 0.0046 
(5£;(^,0) = -0.01478, (5e/ = -0.0056 



(60) 



The band energy is lower for incommensurability along 
the Cu— O bond direction. 

The Cu— O bond direction incommensurability is lower 
in energy due to the additional Umklapp scattering avail- 
able for q on the Brillouin zone edge rather than inside 
the zone for diagonal q. 

We have shown that chiral coupling of polarons to <P 
spins leads to a ring of incommensurability centered at 
(tt/o, 7r/a) of magnitude {2Tr/a)x. The perturbation to 
the kinetic energy of the delocalized x'^ — band elec- 
trons selects incommensurability along the Cu— O bond 
direction due to Umklapp scattering on the Brillouin zone 
edge. 

For 0.02 < a; < 0.05,^^ La2-xSr:rCu04 is a spin- 
glass. No x^ — band is formed because the plaquette 
polarons do not percolate. The — states triplet 
couple to polaron spins. The spin interactions in the 
spin-glass phase are different from the Hamiltonian in 
equation [44l 

The Cu — cannot delocalize over an infinite po- 
laron swath in our model. We do not know if the x'^ — y^ 
states remain localized on a single Cu or delocalize over 
the finite swath of the polaron. Any delocalization leads 
to an effective ferromagnetic coupling between neighbor- 
ing polaron spins due to the triplet coupling with the 
x^ — y^ spin. In addition, there is an asymmetry in the 
chiral coupling due to orthorhombic crystal symmetry 
arising from the tilt of the CuOg octahedra. The one- 
dimensional incommensurability in the spin-glass phase 
of La2-a:Sra;Cu04^ is not explained in this paper. 
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FIG. 14: Effective polaron orbitals with absorbed into Aig 
on the Cu sites. 



VI. SUPERCONDUCTING PAIRING 

Coulomb scattering of — y^ band electrons with chi- 
ral plaquette polarons leads to anisotropic Cooper pair 
repulsion. The maximum energy difference between a 
chiral polaron and its time-reversed partner is analogous 
to the Debye energy in BCS superconductors. As dis- 
cussed in section fVB I this energy difference can be non- 
zero and on the order of Jm- 

We simplify the polaron wavefunctions by absorbing 
the Pa- orbitals into the Aig orbitals on the Cu sites for 
the pairing and transport calculations. This is shown in 
figure UM 

The direct and exchange Coulomb terms coupling the 
x^ — y^ band with a polaron is. 



Hu = U d^B.a'A,aZRadRa'. 



Hk = -K rfflo-'4a2fl<T'C^ 

R,,(T.cr' 



(61) 



(62) 



where U > K > 0. creates an x"^ — electron at R 



and creates an Aig electron at R. 



y 



The x — y band state with momentum k is, 



4. = ^-^E^"'-''4., 



(63) 



R 



where N is the total number of Cu sites. A polaron state 
of spin s is given by. 



4 = Y "R-^R« 



(64) 



R 
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FIG. 15: Second-order Cooper pairing processes with po- 
larons. The Coulomb couphng is represented by a four-point 
vertex for simplicity. The solid lines are the band electrons 
and the dashed lines are the polaron states. / is the initial and 
final polaron orbital state. /' is the intermediate state, (a) 
and (b) are the Coulomb couplings with no spin exchange, (c) 
and (d) are spin exchange couplings. In (c) and (d), the final 
electrons are interchanged compared to (a) and (b) yielding 
an extra minus sign in the pairing matrix element. 



where the coefficients determine the type of polaron 
in figure fHl The matrix elements for direct and exchange 
Coulomb scattering of an x'^ — band electron with a 
polaron electron are, 

{k'a',z^,\Hu\ka,z:) = Af(k' - k),5,,,5,,, , (65) 



M(q)=^a-ai,e-"i^. 



(66) 
(67) 



R 



The Cooper pairing matrix element for scattering aris- 
ing from figure fTSl a) with initial polaron orbital state / t 
is, 



((k't,-k'i),/T |i^'l(kT,-ki),/T> = 



iV2 



Ef — En Ei — En 



(68) 



where q = k' — k, Ei is the total energy of the initial 
state, Ef is the final state energy, and i?„ is the interme- 
diate state energy. The factor U — K comes from the left 



vertex where both direct and exchange Coulomb matrix 
elements appear while the factor U arises from the right 
vertex where there is no exchange. Ei, En, and Ef are. 



E, ^ 2€k + E^, 

En = f'k + efc + -B/'T: 

Ef — 2ek' + Eji . 



(69) 
(70) 
(71) 



Substituting equations [69l [TOl and ITT] into [68l leads to. 



((k' t,-k'i),/T |ff'|(kT,-ki),/T) = 



U{U-K) 
7V2 



|M(q)P5] 



■jf (e/c' - Cfc) 



{Er^ - EiT;) ^^^^ 



The pairing matrix element is the probability weighted 
average of equation [72] over all initial polaron states /. 
States with Ej>cr — Ej^ > dominate the average because 
the polaron is primarily in its lowest energy state. Since 
U > and U — K > 0, there is an attractive pairing 
with momentum anisotropy determined by |M(q)p for 
pairs close to the Fermi level. The coupling is attractive 
for Cooper pairs near the Fermi level within the polaron 
energy separation, |efc| and \ek>\ < {Eji^y ^ Eja-)- There 
are similar terms for / J, and figure [TST b) . 

For the spin exchange processes shown in figures 115( c) 
and (d), the matrix element picks up an extra minus sign 
compared to figures [TSTa') and (b) due to the exchange 
of k' t and — k' j in the final state. 



(-) 



((k'T,-k'i),/^|i/"|(kT,-ki),/a> 
^^'|M(k' + k)P^. ^ 



N 



J, (ffc' - Cfe) 



A2 



(73) 



where A = Eji^ — Ej^. This leads to an overall repulsive 
interaction for pairs near the Fermi level. 

The direct on-site isotropic Coulomb repulsion between 
Cooper pairs must be included in the net pairing inter- 
action in addition to the second-order pairings in figure 

m 

For hole-doped cuprates, one chiral polaron state is 
singly occupied and the remaining polaron states shown 
in figure [3] are doubly occupied. The intermediate po- 
laron state /' must be a chiral state, Px' ± iPy'- For 
electron-doped cuprates, the initial state / is a chiral po- 
laron. The pairing is identical because the matrix ele- 
ment between / and I' is mod-squared. We focus on the 
hole-doped derivation here. 

The momentum dependence of the pairing is deter- 
mined by |M(q)p where q = k' — k for non-spin ex- 
change attraction and q = k' -I- k for spin exchange re- 
pulsion. |M(q)p is the same for / = 5* and Dxy From 
figure [Ml and equation \67\ 



|M(q)| 



2 

S.D 



• 2 1 . 

sm - [qx + qy) 



sin'^i-qx + qy), (74) 



where we have dropped all constants and wavefunction 
normalizations in the expression because they do not af- 
fect the momentum dependence. For / = Fx' ± iPy' and 
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TABLE IV: The values of |A/(i7)|^ at special vectors in the 
Brillouin zone. For non-spin exchange couplings shown in fig- 
ure [I5t^ a) and (b), the last column does not apply because 
the intermediate cannot be identical to the initial polaron 
state, r 7^ /. M{q) is defined in equations [TH [75l and [76l 
q = k' =p k for non-spin exchange and spin exchange respec- 
tively. 
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cos - {q^ + Qy ) 



cos - {-Qx + Qy) 



(75) 

For spin exchange repulsion, there is an additional pos- 
sibility where I' = I = P^' ± iPy'. \M (q)|++ is. 



|M(q)|^ 



cos - {qa: + qy) 



-cos-(- 



qx + qy) 



(76) 

Table HVl shows the value of |M(q)p for various vectors 
q. There are several pairing possibilities. 

First, consider non-spin exchange attractive pairing 
through figures fTSTa) and (b) added to direct isotropic 
repulsion. In this case, q = k' — k. |M(q)|| ^ leads to a 
net attraction for q = (tt, 0). An s-wave gap cannot occur 
because the isotropic repulsion dominates. A gap of any 
other symmetry has a net repulsion for q = (0, 0) and 

cannot lead to superconductivity. |Af(q)|^ also repels 

at q = (0,0) and does not superconduct. M++ does not 
occur for non-spin exchange. Thus, there is no supercon- 
ductivity through non-spin exchange polaron coupling. 

For spin exchange repulsion in figures llSt c) and (d) 
added to direct isotropic repulsion, Ms,d has the largest 

repulsion for q = k' + k = (tt, 0), M-| for q = (tt, tt), 

and Af++ for q = (0,0). Af^^ cannot be correct for 
the cuprates because q = when k' = — k. This leads 
to a gap that changes sign for k — k and triplet spin 
pairing. For Ms^d, there is no favorable Fermi surface 
nesting near k' -I- k = (tt, 0). 

has the time-reversed chiral polaron as the in- 
termediate state and has a large repulsion for k' -|- k = 
(tt, tt). For k sa (tt, 0) and k' k, (0, tt), there is a large 
repulsion between pairs close to the Fermi surface. Pair- 
ing through spin exchange coupling with time-reversed 
polarons leads to d-wave superconductivity. 

As discussed in section |VB] on the neutron spin incom- 
mensurability, the energy splitting between time-reversed 
chiral polarons is on the scale of Jdd ~ 0.1 eV. It is this 
scale that is analogous to the phonon Debye energy in 
conventional superconductors. 

The energy separation between two time-reversed chi- 
ral polarons, \Px' +iPy', a) and \Pxi —iPy', —a), is largest 



for low dopings where there are more undoped (P spins. 
The energy splitting decreases with increasing doping 
leading to a reduction in the energy range surrounding 
the Fermi surface involved in Cooper pairing. Concur- 
rently, the number of polarons that can induce Cooper 
pairing increases. The number of band — electrons 
available to form the superconducting ground state also 
increases with increasing doping. Finally, as the number 
of polarons increases, they crowd into each other, even- 
tually creating polarons localized on single Cu sites, as 
found in UB3LYP calculations^. The single Cu state is 
not chiral and leads to isotropic Cooper pair repulsion. 
The competition between all of these factors leads to the 
superconducting phase diagram in cuprates. 



VII. NORMAL STATE TRANSPORT 



A. Resistivity 



The temperature dependence of the resistivity arising 



from 



y band electrons Coulomb scattering from 



polarons is shown to depend on the density of polaron 
energy separations, /9(A). If the density is constant, then 
the resistivity is linear in T. 

The scattering rate of an x^—y^ electron initially in the 
state kcr to Coulomb scatter from a polaron to another 
band electron and polaron is the sum, 



1 



27r 



J2 \{k'a',I's'\H'\ka,Is)\'P{Is) 



Is, I'. 



X (5(ek'<T' + Eps' - Ekcr - Els), 



(77) 



where H' — Hjj + Hk is the total Coulomb interaction 
of band electrons with polarons. Hjj and Hk are defined 
in equations 1611 and 1621 Eis and Eiis' are the initial and 
final polaron energies respectively. P{Is) is the proba- 
bility the polaron is initially in the state Is. The delta 
function enforces total energy conservation, 



Cko- + Ejs — Sk' 



Er 



(78) 



The matrix element in equation [77] is proportional to 
|Af(k' - k)|2 defined in equations El [73 and [75] To ob- 
tain the qualitative temperature dependence, we simplify 
equation [77] by assuming the matrix element is constant 
for all momenta and there are only two polaron states 
at any doped center. This does not alter the tempera- 
ture dependence, but smears out the anisotropy of the 
scattering rate. The M{k' — k) have different momen- 
tum dependencies for each polaron symmetry. Since the 
scattering rate is a sum over all scattering channels, the 
approximation of isotropic scattering is likely to be good. 

There is a polaron Coulomb scattering process we have 
neglected. This is second-order hopping of a k cc^ — j/^ 
electron into the unoccupied polaron state that then hops 
off into a k' — band state. Since the polaron is 
comprised of orbitals with Ai symmetry, the hopping is 
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zero for momenta along the zone diagonal and largest at 
(±7r, 0) and (0, ±7r). This term is approximately temper- 
ature independent and of "d-wave" symmetry. The tem- 
perature dependence of this term arises from the displace- 
ment of the apical O due to the addition of an electron 
into the polaron. Since the initial and final states in this 
second-order process have one hole on the polaron, the 
temperature dependence is expected to be weak. We ne- 
glect this additive temperature independent anisotropic 
constant term in the remainder of this section. 

The neglected term justifies the addition of a d-wave 
lifetime broadening of the band electrons for the compu- 
tation of the change in the band energy from the neutron 
spin incommensurability in section fV CI 

Consider an — band electron with energy ei and 
a polaron with energy Ei scattering to an electron with 
energy £2 and polaron with energy Ef. ei and e2 are mea- 
sured relative to the Fermi level, = 0. Conservation of 
energy is, 



€1=62 + A, A EE £: 



/ 



(79) 



The thermal occupations of a doped site with two po- 
laron states of energy Ei and Ef is given by the Fermi- 
Dirac functions f{—A) and /(A) defined in equation [24] 
and depends only on the energy difference A. This is 
identical to the results in equation 1231 

The total lifetime of a band state with energy ei is 
given by, 



1 



J de2 J dA p{A)'^\H' 



X <5(e2 - ei - A)/(-A) [1 - /(e2)] , (80) 



where p( A) is the density of polaron sites with an energy 
splitting A and is taken to be a constant. The p(A) 
used here is different from the spin and chirality flipped 
density in section fVAl p(A) in equation [501 includes the 
energy splittings of the S, D^y^ and the bonding combi- 
nations of the Px' and Py' polaron states. 
Integrating over the final electron energy €2 , 

1 = ^|i/?iV(0) 
T(ei) n 

X 1 dA p(A)/(-A) [1 - /(6i - A)] , (81) 

where iV(0) is the band density of states per spin. 

For ei at the Fermi level, ei = 0, the integrand becomes 
p(A)/(A)[l — /(A)] leading to a temperature dependence 
T^^'^^ if p(A) ~ A^. For a uniform distribution of polaron 
energy splittings, p(A) — po for |A| < A^ax and zero 
otherwise, the integral in equation [SD is. 



27r, 



r(ei) 



H'\'NiO)po 



T 



(1- 6-/5^1) 



In 



„/3A„ 



-0ei 



-^1) + 1 



(82) 



The energy scale for Amax is determined by the scale 
for the hopping matrix elements and is several tenths 
of an eV. This is much larger than the temperature, 
f3Amax S> 1. For |ei| ^ Amax, wc may take the limit 
Amax +00 leading to. 



ei 



1 ^^|H'|2^(0)p 

r(ei) n V 1 — e 



(83) 



The tb that appears in the Boltzmann transport equa- 
tion is related to the lincwidth and satisfies, 



Ts(ei) h r(ei) l-/(ei) 



27r 



1 



TBie) 



= —\H'\'NiO)po ecoth-/3e. 



At the Fermi level, e = 0, the scattering rate is. 



1 



2tt. 



H'\'N{Q)po{2T), 



(84) 



(85) 



(86) 



leading to a linear resistivity. 

If p(A) ~ A for |A| < Amin and is constant for 
Amin ^ A < Amax, then the resistivity is ~ at 
low temperature and crosses over to ~ T at high tem- 
perature. Different forms for p(A) lead to different 
temperature dependencies for the resistivity as seen in 
La2-:,Sr:,Cu04 from 0.05 < a; < 0.35.8° 

The optical spectrum is composed of a Drude scatter- 
ing term for finite e in equations [82] and [85] plus polaron to 
polaron scattering due to photon absorption. The Drude 
scattering rate becomes proportional to e for /3e ^ 1 with 
the crossover from T to e occurring when /3e « 2. When 
the energy is larger than the largest polaron splitting, 
e S> Amax the scattering rate saturates and becomes 
proportional to Amax- 

Direct optical absorption from S and D^y polarons 
to Px' ± iPy' can occur for light polarized in the Cu02 
planes. This leads to the excess mid-IR absorption;2i 

Chiral polarons can lead to the logarithmic resistivity 
upturn at low temperatures due to a Kondo effect. As 
La2_2;Sr2;Cu04 is doped, polaron islands are formed com- 
posed of exactly one polaron with no adjacent polarons. 
Island polarons are surrounded by spins. They have 
two possible chiral ground states. Antiferromagnetic spin 
flip scattering with the — band occurs in second- 
order where the — y^ electron hops onto the polaron 
and either an up or down spin polaron hops back to the 

— y^ band. This leads to a Kondo effect. 

The momentum dependence of the antiferromagnetic 
Kondo spin exchange is largest for momenta near (±7r, 0) 
and (0, ±7r) since the polaron is comprised of Ai orbitals. 
The coupling is zero for momenta along the diagonal. 
The resistivity in the plane is dominated by electrons 
with momentum near the diagonals and transport out 
of the plane by momenta near (±7r, 0) and (0, ±7r). The 
resistivity upturn due to Kondo scattering appears at a 
higher temperature for out of plane transport. 



21 



0.62S 




O {f.{l2 {f.&4 0.06 OffS OJ OJ2 SJ4 OJ6 £>.Ji 

Doping 

FIG. 16: Number of polaron islands per Cu as a function of 
doping for a 2D Cu lattice. An island is defined as a four-site 
polaron having only undoped d? spins as neighbors. 



There is no Kondo effect from Coulomb scattering with 
polarons from Hu + Hk defined in equations [H] and [M] 
because the spin coupling is ferromagnetic. 

Figure [1^] shows the number of islands per unit cell 
as a function of doping for a 2D lattice with doped po- 
larons. The Kondo resistivity is expected to disappear as 
the number of islands goes to zero. From table IIIIl the 
polarons percolate in 2D at x « 0.15. This is approxi- 
mately the doping where the insulator to metal transition 
occurs for La2_2;Sr3;Cu04.®^ 

The above suggestion has several caveats that could in- 
validate the conclusion. First, we do not know how many 
islands exist with a small energy separation between the 
chiral polaron time-reversed states. 

Second, the g-factor for these islands needs to be eval- 
uated to determine the energy splitting between the "up" 
and "down" states. The Kondo effect is suppressed for 
temperatures less than the splitting energy. 60 Tesla 
pulsed magnetic fields were used to suppress the super- 
conductivity in order to access the low temperature nor- 
mal state 1^ 

The g-factor of an electron spin is two, = 2. A 60 
Tesla field splits the up and down electron spin energies 
by 80K. For g = 0.1, the splitting is 4K. The Kondo effect 
is active for T > 4K. 

An estimate of the polaron g-factor, gp, is made in sec- 
tion IVII Bl by fitting to the temperature dependent Hall 
effect for La2_a;Sr2,Cu04 at a; = 0.10. It is argued that 
gp ~ 0.1. Also, for island polarons, figure [TOl shows the 
energy difference between a polaron and time-reversed 
polaron is zero for large Jet- This order of magnitude 
reduction in gp allows the Kondo effect to remain active 
at low temperatures. 

Third, the coupling to the band electrons needs to be 
evaluated to determine if the logarithmic resistivity is of 
the right magnitude. 



B. Hall Effect 

Skew-scattering has been proposed^^i to explain the 
temperature dependence of the Hall effect although the 
physical nature of the excitations causing the skew- 
scattering was unclear. In this section, we derive and 
estimate the magnitude and temperature dependence of 
the skew-scattering of — band electrons from chi- 
ral plaquette polarons. The skew-scattering contribution 
to the Hall effect is an extra term that is added to the 
ordinary Hall effect. 

Skew-scattering2^i21iS^iSi2i is a left-right scatter- 
ing asymmetry occurring when the scattering rate from 
k — > k' is not equal to the scattering rate from k' k, 
^(k k') w(k' k). A combination of time-reversal 
and inversion symmetry leads to no skew-scattering be- 
cause w(k — > k') = w(— k ~* — k') due to inversion and 
k — > — k') = ti;(k' — !■ k) from time-reversal invari- 
ance. An applied magnetic field breaks time-reversal in- 
variance by making the number of polarons of "up" and 
"down" chirality different. 

Skew-scattering first appears in third order for the 
scattering Hamilitonian as can be seen in the simple ex- 
ample below where wc ignore the polarons and consider 
a band scattering matrix element that is complex. 

Let H' be the electron scattering Hamiltonian, 



(87) 



k,k'o- 



where we ignore any spin dependence in the matrix el- 
ement Vkk'- Since H' is Hermitean, 14;/^ = ^k" The 
scattering T matrix is given by 



T{E) = V + VG{E)V + ■ ■ ■ , 
and the scattering rate is, 

wik^k') = ^\{k'\Tie,)\k)\'5ie,,-e,). 

Expanding the T matrix to second-order, 

^k'k"^k"k 



(89) 



(k'|r(6,.)|k) = Fk-k + E 



t;' cfc - cfc' + i5 



0{V^). (90) 



Substituting into equation [55] and neglecting terms of 



w{k k') 



2tt 



iMc'kP 



E 



2Rc(Fk.kVkk"V^k"k') 



-27r^ Im(Vk'kVkk"Vk"k')<5(efc - ffe") 

k" 



6{ek-ek'). (91) 



As a check, equation[ni]is invariant under any redefinition 
of the k states, Vk'k ^ e'(^i"-ek)y^,|^ 

If V is real, then there is no skew-scattering since 
^(k ^ k') = w{k'^k). If F is complex, then in- 
terchanging k and k' changes the sign of the third 
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term. This is the lowest order skew-scattering term. 
By interchanging k' <-!■ k" in the third term of the 
equation, the sum over k' of the skew term satisfies 

Sfc' WskewO^ k') = 0. 

The applied magnetic field causes skew-scattering by 
creating left-right asymmetries. The relevant question 
is not if there is any skew scattering, but whether the 
scattering is large enough to account for the experimental 
Hall effect and its temperature dependence. 

We derive the skew-scattering from the Coulomb 
Hamiltonian Hu in eauation lGll Including the exchange 
term Hk in equation [B^l does not change the results be- 
low because it leads to an average over U and U — K. 
Thus, we focus on Hu and ignore the electron spin. 

Consider a single polaron. In the hole-doped cuprates, 
one of the chiral polaron states, = Px' ± iPy' is ini- 
tially unoccupied and all the remaining polaron states 
shown in figure [TJ] are occupied. The matrix elements 
for scattering an — k electron and polaron electron 



The T matrix element is defined as. 



(100) 



where the hole resides in (the electron is in P^) 
Expanding T^j, to order U'^, 



k'k 



k'k 



(101) 



where the sum is over all intermediate polaron states /. 



^k'k defined as. 



yd) 



k'P±|i?r/|kP=' 



(102) 



and is independent of the polaron chirality and real from 
equation [Ml Thus, V^l = V,^^^'- ^^^^ is. 



to a chiral polaron state with a change in the 
band electron momentum to k' is. 



y 



T/(2)(±) 



k'P^\Hu\pI){pI\Hu\kP^) 

ek-ep-{Ei-E±) + i5 ' 



(103) 



{k' P^\Hu\kS) 



k:P^\Hu\kDxy) 



(k'P±|iJr/|kP±) 



U 
2N 



U 
2N 



{-i) (sin (7a;/ ± isinqyi) 



where the sum is over all intermediate — momenta 
(— i) {sinqx' T isinqyi) , (92) ^± g^j^j jtj^ ^j^g chiral polaron and the intermedi- 

ate polaron energies respectively, ep is the — y^ band 
energy for momentum p. 
Since V^^'> is real, \T^,^\^ is, 

iT^k^kP = {vSiy + V^iH E 2RcVSl^^\l) + 0(C/3). 

(104) 

The skew-scattering contribution is antisymmetric under 
interchange, k' ^ k. The first term in equation 11041 is 
symmetric and so is V^f^^^\l) for / = P^. The lowest 
order skew-scattering term is. 



U 
2N 



(93) 



{cosQx' + cosqy') , (94) 



(j^^ cos ^q^ cos ^qy, (95) 



(k'P^lHulkP^) 



where q = k' — k and 



U 
2N 



(cos qx 



cos qy' 
. 1 



(96) 



l^k^kl^ = 4'^2Re 



vSl^^\s) + vSl^^\Dxy)]. (105) 



(-) ( ^ j sinig^sin^gy, (97) 



In appendix |B1 it is shown the antisymmetric terms in 
V(^^(^^{S,Dxy) are, 



i-Qx + qy) ■ 



(98) 
(99) 



(106) 

vrPi [ek-{ED-E^)] Ak',k, 



The chiral polaron to chiral polaron scattering terms 
in equations [15] and [§7] are real and lead to no skew- 
scattering, as shown below. Second-order mixing with S 
and Dxy polaron symmetries leads to skew-scattering. 

The Coulomb repulsion, U, in the above matrix ele- 
ments is positive, U > 0. The matrix elements corre- 
sponding to equations [^^H^ for hole polarons are ob- 
tained by the substitution U ^ —U. We use the expres- 
sions in equations [5^1-1571 and remember to take U < 
for the hole-doped materials and C/ > for the electron- 
doped materials. 



El 

m 



(107) 



where the antisymmetric function, Akk' = ~^k'k, is, 

^k',k = sin fca;' sinky, — sin fcj,/ sinfc^/. (108) 

Equations [Mj and [Mj define the x' and y' components of 
k and k'. The function Pi is the sum over the Brillouin 
zone, 

Pi(w) = E ^ (cospa: + cospy) S{uj - ep). (109) 
p 
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The Dxy contribution is negative the S contribution with 
Eo substituted for Es- 

The skew-scattering term is, 



2tt 



1 



cos -q.j; cos -Qy 



— 1 (- 

V 2 

{Fi [eu - {Es - E^)] - 
Fi[tk-{ED-E^)]]Ai,,^, 



1 



(110) 



where q = k' — k. 

Equation lllOl is the skew-scattering for a single polaron 
with polaron orbital energies E^, Es, and Ed- These 
energies have probability distributions p'^{E'^), ps{Es), 
and pd{Ed)- The mean value of equation [110] is, 



dE^dEsdEDP^{E^)ps{Es)pD{ED) 

X \T^,^\l{E^,Es, Ed). {Ill) 



If polarons with S and Dxy symmetry have identical en- 
ergy distributions, then the ± components of the skew- 
scattering are always zero. 

The range of energies that E'^ and E~ span is ~ Jdd 
and is much smaller than the hopping matrix element 
scale of Es and Ed ~ tuop- Therefore, we set E^ = in 
Fi and evaluate the mean over Es and Ed, 



<.eJk - k') = -{\T^.^\l)5{eu: _ eu). (116) 



The total skew-scattering is, 

WskeA^ - k') = N.p{n- - n+)7«+,„(k ^ k'), (117) 

where Np is the total number of polarons and {nT — n^) 
is the mean number of P~ polarons minus the number of 
polarons. This difference is non-zero in a magnetic 
field and is proportional to the field B. 

To evaluate (n^ — W^), consider a polaron with energy 
difference between P" and P+ of A, E+ - E- = A. The 
field changes the energies to E^ = E^ ^ gpPBB, where 
gp is the polaron g-factor and ps is the Bohr magneton. 
The occupations become. 



n ^ fi-A) + 2gpPBB{-^ 



df 



B + 0{B^), (118) 



n+ = /(A) - 2gpPBB 



dl 
V de 



B + 0{B^), (119) 



where /(e) is the Fermi-Dirac function defined in equa- 
tion [24l Taking the probability density of A to be uni- 
form over the range Ami„ < |A| < lS.max makes the 
density po = l/2(Amaa: — Amin). The mean polaron dif- 
ference is. 



{Fi[ef~Es]) = / dojp''{oj)Fi{€f-oj), 
Jo 



AS ' 

max 



(112) 



where 



1 



^{cosp^+ cos py), (113) 



and pS = 1/A;^ 



the interval < A < A£„^ and 



zero outside. There is a similar expression for the mean 
over Ed- 

If £/ — A^£, is less than the bottom of the x^ — y"^ band, 
then A"^ = = A where A is the sum over occupied 
states, 



{cOSPx + COS py). 



(114) 



In general, — 1 < A < 1. Cuprate band structures are 
occupied at (0, 0) and unoccupied at (tt, tt) leading to 
< A < 1. For La2-a;Sr:rCu04 at a; = 0.10, A = 0.19. 
Substituting the mean of Fi into equation II 101 



TT^ \ U [ — ] COS -Qx COS -Oy 

I \2) 2 2^ 



X A 



1 



1 



A^^ 



A^ 



k'k, 



(115) 



(n - 71+) = AgpPBB(2po) 



dA{-^), (120) 



{n - n+) = AgpPB 



/ (^maa; )] 



B. (121) 



All of the temperature dependence of the Hall effect 
is contained in equation 11211 Amax is on the order of 
Jdd ^0.1 eV. Ajnin and Amax are larger for small doping 
since there are more undoped d^ spins to split the energy 
between time-reversed polarons. The Amax in equation 
11211 is smaller than the Amax in equation [52] since the 
latter includes the splittings of the S, D^y, and bonding 
combinations of P^' ± iPy'- From equation I121i the dif- 
ference (n^ — n+) is zero for T <^ A,„j„ and rises to a 
maximum for some T between Amin and Amax- When 
T Amax, the difference decreases to zero as l/T. The 
temperature dependence from skew-scattering is added 
to the ordinary band contribution to the Hall effect. The 
ordinary term is hole- like and positive for the cuprates. 
The coefficient that multiplies equation 11211 in p^y'^'^ can 
be either positive or negative depending on the details 
of the polaron energy distributions and the sign of the 
Coulomb repulsion U . 

The observed temperature and doping dependence of 
La2_a;Sra;Cu04^ is consistent with chiral plaquette po- 
laron skew-scattering as shown in figure [T7] We fit the 
experimental data to the expression. 



Rh = Ro [/(A n 



(122) 
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FIG. 17: Fit to Hall effect dataS' for La2-^Sr;rCu04 using 
the expression in eauation ll22l For x = 0.05, the parameters 
are A™i„ = 12.6 K, A^a^r = 791.5 K, and 7?o = 15.4 x 10"^ 
cm^/C. For x = 0.10, A„i„ = 16.7 K, A„„^ = 454.1 K, and 
i?o = 8.82. For x = 0.15, A™i„ = 0.0 K, A™a^ = 387.9 K, 
and i?o = 3.81. For x = 0.25, A™„ = 11.6 K, Amax = 106.3 
K, and Ro = 1.05. 



where /(e) is the Fermi-Dirac function defined in equa- 
tion [21 

There is sufficient structure in the temperature depen- 
dence derived above to account for the various tempera- 
ture behaviors of the electron-doped materials, tooi^ 

The additive contribution to the Hall resistivity arising 
from skew-scattering is derived in appendix [Bl and is. 




To determine if the magnitude of the chiral polaron 
skew-scattering is large enough to be compatible with 
experimental data, we computed RhC = p'^^y^/B for 
La2-2;Sra;Cu04 using the band structure in equation 1581 
at X = 0.10 and estimates of various parameters. 



We take m to be one electron mass and n = Ax/^ceii 
to be the number of charge carriers per volume (since 
each polaron adds four — band electrons in our 
percolating model). The polaron g- factor is taken to be 
one, gp = 1. A^ax = 454.1 K and Amm = 16.7 K from 
the fit in figure flTl Finally, we choose A^^^ = 1.0 eV 
and ^max = 1-5 eV for the energy distribution widths of 
the S and D^y polaron energies. 

The reason we chose A^^^ > A^^^^ is because the 
energy separation of the D^y states from the chiral 
polarons is larger than S in table HH This does not prove 
^max > ^fnax always the correct choice, but it is sug- 
gestive of the expected energy splitting. Also, our values 
for the average skew-scattering T matrix in equation [115] 
assume a uniform energy density and p^ and this may 
not be correct. The idea here is to get a reasonable order 
of magnitude estimate for the size of the skew-scattering. 

The final parameter to choose is the Coulomb repul- 
sion matrix element energy U. [/ ~ 1.0 eV is the cor- 
rect energy scale because U — {l/nceu)4:ne'^ /qj-, where 
qd ~ 47re^7V(0) is the Debye screening length arising 
from the many-body response of the band electrons to a 
perturbing potential. A^(0) is the band density of states. 
A typical metallic density of states is one state per eV 
per unit cell leading to C/ ~ 1.0 eV. Since the density of 
charge carriers is reduced to Ax = 0.40 per unit cell, the 
band electrons are less effective at screening the Coulomb 
repulsion. This increases C/ by a factor of 1/0.4 = 2.5. 
We take U = 2.5 eV. 

U is negative for the hole-doped cuprates because the 
band structure arises from fully occupied polaron states 
and a hole in a chiral polaron state amounts to subtract- 
ing a Coulomb matrix element. This was discussed previ- 
ously following equations [92l-l97l For the electron-doped 
systems, U is positive. Since appears in equation 
11151 the sign difference in the temperature dependence 
of the Hall effect between hole-doped and electron-doped 
cuprates is obtained by assuming all the other parame- 
ters discussed above do not change. We take U = —2.5 
eV for La2-a;Sra;Cu04. 

Calculating p^y'^^ in equation 11231 we obtain Rhc — 
5.39 X lO-^cm^/C at lOOK or a charge density of 0.116 
carriers per unit cell with the chosen parameters. The 
experimental value-- for x = 0.10 at 100 K is approxi- 
mately 4 X 10^'^cm'^/C as can be seen in figure [T71 Our 
result is a factor of 1.35 too large. 

The largest errors in our calculation come from the 
estimate of the Coulomb repulsion U, the effective 
(n/m)eff, and the value of polaron g-factor, gp = 1. We 
used one electron mass for m and chose n to be 4a; = 0.4 
carriers per cell to obtain (n/m)eff. This is likely an 
overestimate of {n/m)ef / since the carriers must traverse 
a percolating swath through the crystal and the scatter- 
ing rate 1 /r for the ordinary conductivity is not expected 
to be large. Since {n/m)~^f and appear in p^y^ small 
changes lead to a large enhancement of the computed 
Pxy We conclude the polaron g-factor is overestimated 
and may be an order of magnitude smaller, leading to 
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(7p ~ 0.1. A reduced gp allows the Kondo resistivity dis- 
cussed in section IVII Al to remain active in the pulsed 60 
Tesla magnetic fields^ used to suppress superconductiv- 
ity and measure the low temperature resistivity. 

The analysis of this section can be done for the real 
polaron states Px' and Py' instead of the chiral states. 
In this case, the magnetic field creates the Van-Vleck 
paramagnetic states Px' + iBPy> and Py' + iBPx'. The 
B in the equations is a constant times the applied field. 
The skew-scattering contribution from a polaron where 
the energy of Py' is greater than Px' by A exactly cancels 
the contribution from a polaron where Px' is greater than 
Py' by A. Thus, real polaron states cannot lead to skew 
scattering. 

For a planar magnetic field, the g-factor is substan- 
tially smaller because the angular momentum of the chi- 
ral polarons is normal to the plane. Therefore, for planar 
magnetic fields, the Hall effect should be temperature in- 
dependent, as observed. 



VIII. ARPES 
A. Pseudogap 

The undoped spins have antiferromagnetic fluctua- 
tions with wavevector q « (tt, tt) that couple an — 
band state with momentum k to k ± q. In this section, 
we show that this coupling leads to a pseudogap at the 
Fermi energy with magnitude proportional to the square 
of the couphng. 

At fixed temperature, the coupling is larger for low 
dopings because there are more undoped (P spins. The 
couphng decreases as the doping increases. This leads to 
a larger pseudogap for low dopings, as observed i^ii^^i^i^i 

At fixed doping, the strength of the antiferromagnetic 
(fi fluctuations decreases with increasing temperature. 
Thus, the coupling to the x"^ — band should decrease 
with increasing temperature. This leads to the closing of 
the pseudogap with temperature. 

The momentum dependence of the pseudogap is deter- 
mined by the shape of the Fermi surface. If the state 
with momentum k ± q is close to the Fermi surface, then 
mixing with the k state leads to a reduction in the spec- 
tral function and a pseudogap. Such nesting occurs for 
states near (±7r, 0) and (0, ±7r). Along the Brillouin zone 
diagonals, there is no nesting. This leads to a pseudogap 
that is zero along the diagonal and increases as one moves 
toward (7r,0). The model does not guarantee zero pseu- 
dogap along the diagonal although the magnitude of the 
gap is expected to be smaU. 

ARPES measures the spectral function ^(k, cj) = 
(— l/7r)Im G(k, w). Evaluating A(k, w) for k vectors on 
the Fermi surface, Ck = £/, and at the Fermi energy, 
Lj = €f, determines the density of states suppression at 
the Fermi level and its temperature dependence. The 



relative change in the spectral function is, 

.(k...V)= -^"--"'.-'f"' . ,124, 

where Ao(k, lS) is the V = Q spectral function and A(k, uj) 
is the spectral function including the (fi interaction V . 
The pseudogap is zero when V ~ Q &i high temperature. 
Using equation ITZl 



^^y2, Im C;g(k,^ ) [Go(k + q, + Go(k - q, w)] 



Im Go(k, Lo) 



From equation I51[ 



Go(k,eO 



iFt 



(125) 



(126) 



The relative spectral function change at the Fermi level, 
uj = e/, is. 



V 



k+q 



Ffc-g 



{ef-ek^qr+Tl_^\ 



(127) 



r is always negative, r < 0, and r —>■ as V ^ 
at high temperatures. The density of states decreases 
at the Fermi energy leading to a temperature dependent 
pseudogap that closes as the temperature increases. The 
linewidths F^ and Tk±q are of the order ~ 0.01 eV due to 
the broadening of a k state arising from the non-uniform 
percolating path in underdoped systems. The coupling 
of the dp spins to the x'^ — band must be less than 
or of the order of Jdd ~ 0.1 eV. The relative density 
of states change is on the order of r ~ (V/Tk)'^- Thus, 
the change in the ARPES spectra arising from coupling 
between the undoped spins and x^ — y^ is large enough 
to qualitatively account for the observations. 



B. 2212 Bi-Layer Splitting 

There is an ~ 0.1 eV splitting between the bonding and 
antibonding x'^ — y^ bands in bilayer Bi2Sr2CaCu208+5 
for k vectors near (tt, 0) arising from coupling through 
out-of-plane orbitals. The splitting along the diagonal is 
very small and not experimentally resolvable. The (tt, 0) 
splitting has been observed^^ by ARPES for overdoped 
systems. No splitting has been observed for underdoped 
systems. In addition, the bilayer splitting disappears as 
the temperature is increased with the temperature sepa- 
rating the coherent (splitting) and incoherent (no split- 
ting) regimes increasing with further overdopingi^ 

XAFS^i^ and flrst principles calculations^ find the 
distance between an apical O and planar Cu decreases 
« 0.2 — 0.3 A due to Sr doping in La2-a;Sra;Cu04. Co- 
herent hopping from one Cu02 layer to another requires 
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coupling of — with an out-of-plane orbital. The 

— band state is delocalized over the percolating po- 
laron swath. Hopping onto an out-of-plane polaron or- 
bital induces a local structural deformation. The matrix 
element is reduced due to the large structural change as 
in the small polaron problemi^ The — band elec- 
tronic response to the changed potential of the doubly 
occupied polaron suppresses the matrix element due to 
the orthogonality catastrophe^^ 

The second-order coherent bilayer hopping through an 
out-of-plane polaron orbital is not completely suppressed 
by the orthogonality catastrophe because the initial and 
final polaron states have the same occupation. Bilayer 
splitting of the — y^ band is induced by this coupling. 
The splitting is strongly temperature dependent and its 
magnitude decreases with increasing temperature. 

In section IVII Al a similar second-order hopping lead- 
ing to d-like anisotropic scattering was argued to be 
weakly temperature dependent. The difference between 
the previous intra-layer and current inter-layer process is 
that the matrix element for the latter is very sensitive to 
any out-of-plane atomic displacement leading to a strong 
temperature dependence. 

Qualitatively, for underdoped 2212, the effective inter- 
layer hopping matrix element is small due to the reduced 
number of Cu pairs between layers that are both part 
of the percolating band swath. In this scenario, the bi- 
layer splitting is small and may not have been resolved 
by ARPES at low temperatures. 

As the doping is increased, adjacent Cu pairs between 
layers become available for coherent hopping and bilayer 
splitting becomes resolvable. This splitting decreases 
with increasing temperature in complete analogy to the 
small polaron problem. At high temperatures, the split- 
ting disappears. 

For highly overdoped materials, the splitting is large 
at low temperatures. As the polaron plaquettes become 
more crowded, any atomic distortion from hopping to or 
from the out-of-plane orbitals is reduced. This leads to a 
higher temperature before bilayer splitting becomes too 
small to be resolvable by ARPES. 



IX. NMR 

A potential problem with our chiral plaquette po- 
laron theory of cuprate superconductivity involves the 
NMR experiments by Takigawa et a& on underdoped 
YBa2Cu3 06.63 that shows the same temperature de- 
pendent Knight shift for both the planar Cu and O. 
The most reasonable conclusion from this experiment 
is that only one electronic component is involved in 
the cuprate NMR. Similar conclusions were reached for 
La2-2:Sr2:Cu04i^i^i2S This would appear to contradict 
our assertion of an out-of-plane polaron because one 
would expect the polaron to have different hyperfine cou- 
plings to the planar Cu and O atoms than the — y^ 
band, leading to a total Knight Shift temperature depen- 



dence (sum of — y^ band plus polaron) that is different 
for Cu and O. 

On the other hand, small-tip angle spin-echo double 
resonance (SEDOR) measurements^" find the electron 
mediated nuclear-nuclear couplings between neighboring 
Cu-0 and Cu-Cu are incompatible with one-component 
cuprate theories. In addition, the high-temperature Y 
spin relaxation and Knight Shifts^^ are incompatible with 
one-component cuprate theories. Indeed, recent phe- 
nomenological fits?2, imply there are two contributions 
to the Redfield correlation time. 

Clearly the effect of the out-of-plane orbitals in our 
model cannot be ignored in interpreting the NMR data. 
In our model, the polarons lead to a weakly temperature 
dependent Knight shift of the form in equation 11211 for 
the magnetic susceptibility due to polarons and shown 
in figure [17] for the temperature dependence of the Hall 
effect. The polaron hyperfine coupling to the planar Cu 
and O sites is weaker than the — y^ band because the 
hole character is primarily in the out-of-plane O Pz or- 
bital. This weak temperature dependence is seen in fully 
doped YBa2Cu307 because the x'^ — y^ band electrons 
Knight shift contribution is constant with temperaturei^i 
For underdoped samples, the temperature dependence of 
the polaron shift is masked by the larger temperature 
dependence of the Knight shift from the x^ — y^ band 
electrons. The latter shift arises from the decrease in 
the density of states of the x^ — y"^ band due to the cou- 
pling with the undoped Cu cP spins as shown in equation 
ll27l for the ARPES pseudogap. We intend to develop the 
methodology to estimate the size of these effects and their 
detailed temperature dependence in order to compare the 
results with experimental data. However, at the present 
time the NMR must be considered a possible Achilles 
Heel for our chiral plaquette polaron theory of cuprate 
superconductivity. 



X. CONCLUSIONS 

Ab-initio quantum mechanical calculations find local- 
ized holes in out-of-plane orbitals in contrast to the t-J 
model. Based on the results of calculations with explic- 
ity doped Sr in La2-a;Sra;Cu04, we postulated that chiral 
plaquette polarons are created by doping and that a de- 
localized x^ — y^ Iva band is formed when the plaquette 
polarons percolate through the crystal. 

The D-wave superconductivity, temperature depen- 
dent Hall effect, neutron lj/T scaling, neutron in- 
commensurability, resistivity, the doping value for the 
insulator-metal phase transition, optical absorption, low- 
temperature log(r) resistivity, ARPES pseudogap, and 
the evolution of bi-layer splitting in Bi-2212 are explained 
by chiral plaquette polarons. 

We have not shown ab-initio that chiral plaquette po- 
larons are formed with doping. For La2-a;Sra;Cu04, ab- 
initio evidence for holes in out-of-plane orbitals has been 
demonstrated. Out-of-plane hole orbitals are plausible 
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for cuprates where the doping arises from interstitial O 
atoms. YBa2Cu307, where the Cu— O chain is com- 
pletely full, is most intriguing. We believe there are local 
structural deformations of the apical O atoms favoring 
localized out-of-planc orbitals. 

Detailed calculations of the spectral function for the 

— band measured in ARPES are necessary in order 
to study the evolution of the background and linewidth as 
a function of doping and temperature. For underdoped 
systems, the percolating path is not uniform through the 
crystal and leads to broadening of the k state spectral 
function. Mixing with the Ai orbitals of the polarons 
increases the broadening near (tt, 0). 

A calculation using the Coherent Potential Approxi- 
mation (CPA)^^'^'' may be the way to obtain quantitative 
estimates of the ARPES. CPA calculations would lead to 
an accurate determination of the magnitude of the resis- 
tivity and would allow us to determine the correct value 
of {n/m)eff to use in the results for the magnitude of 
the skew-scattering contribution to the Hall effect. 

Since neutron spin scattering couples to the electron 
spin and transport couples to its charge, we regard our 
results to be strong evidence for the existence of chiral 
plaquette polarons. 
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The first term on the right-hand side can be expanded in 
a power series in 5ef, 



d 



^no(k, ef)+Sef — 



= ntot+SefNo{ef) + OiSe}), (A4) 

where we have used equations [50l and [53l Afo(e/) is the 
unperturbed density of states per spin. The second term 
in equation IA3I becomes. 



V 



(A5) 



Since 6ef ^ 0{V^), the second term is of order Sej ^ 
O(y^) and can be neglected to lowest order. Substituting 
into IA3I leads to, 



5e/iVo(e/) + l/2 / ' dco ^ P(k, w) = 0. 

J —oo 1. 



(A6) 



From equation IA2| IA6I is identical equation [551 

Substituting lAll into equation [54] for the total energy, 
using equation I50[ and keeping terms up to 0(V^^), 

Etot{(l,V) = Ea + ef -ScfNoief) 



-V 



[ ' cj^P(k,tj). (A7) 



The second term proportional to Stf can be eliminated 
by using equation I A6I leading to, 

Etot{ci,V) = EG + V^ r dc^(c^-e/)^P(k,c^). (A8) 

J — OO 1. 



APPENDIX A: PERTURBED BAND ENERGY 

In this appendix, we derive the expressions used in sec- 
tion [VCj to compute the contribution of the — band 
energy for different incommensurate neutron q-vectors. 

Using equations [351 [STJ and [511 



Eg is the unperturbed {V = 0) ground state energy. 
Using [X2I this is the result shown in equation [551 

To compute the change in energy due to incommensu- 
rability, we need to evaluate integrals of the form. 



IiE,E',ef) 



doj 



00 (t^ - EY {uj - E') 



n(k,e) =no(k,e) + uM dwP(k,cj), (Al) 

— 00 

where we have defined, 

P(k, Lo) = P(k, k + q, w) + P(k, k - q, cj). (A2) 

The total number of electrons is ritot = X^k '^C^' ^f~^^^f)- 
From equation lAll 



J2MK^f+Sef) + V^ r ^dtj^P(k,w). 

1. J — CO 1. 



ntot 



(A3) 



J{E,E',ef)= r diu ^^^^^ , 

i-00 {lj - Ef {ij - E') 

where 

E = e-iT, E' = e' - iV and F, F' > 0. 
From equation I57i 

P(k, p, uj)dLj = 



Im I{ek - iTk, ep - iVp, e/). 

TT ' 



(AlO) 



(All) 



(A12) 
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The imaginary part of J{E, E' ,ef) appears in[ 

The integrand of IA9I can be expanded in partial frac- 
tions, 



1 



{uj - Ef [uj - E') 



= 5E- 



1 



1 



5E-^- 



uj-E' uj-E 
1 



{c-eY 



(A13) 



where 5E = E' — E. The second term on the right-hand 
side can be integrated, 



dw- ^ ' 



1 



{cj - Ef £/ - E 



(A14) 



Each individual integral in the first term is infinite. The 
integral is performed by substituting —A for — c» and 
taking the limit A — > +oo. 



Re / duj 

-A 



1 



1 



i^-E' oj-E_ 



{u - e') + r'2 {u - ey + r2 



(A15) 



Im / duj 

'-A 



duj 



1 



UJ- E' oj-E 



- ef + r2 {lo- e'f + r'2 



(A16) 



The real part can be solved, 



(-A-e')/r' 2;2 + 1 J(-A-e)/r 



■In 



(-A-er/r2 + l 
(-A - e')Vr'2 + 1 



(A17) 
(A18) 



Taking the limit A 
1 



Re / duj 

-oo 



Lo-E' UJ-E 



-oo, the real part becomes, 



■In 



r2 (.t2 + 1) 

(A19) 



The imaginary part integrates to, 
1 1 



Im / duo 



UJ- E' UJ- E 



= tan X 

- tan"^ x', (A20) 



where x — (e/ — e)/r and x' — (e/ — Substituting 
into equation I A91 



r'2 x'^ + 1 



IiE,E\ef) = SE-U-\n 



r2 (a;2 + 1) 
i (tan^^ X — tan^^ x') } 
1 



dE~^ ■ 



ef-E' 



(A21) 



where E and E' are defined in equation lAllI and SE 
E' -E. 



Writing uj-ef = {uj - E') + {E' - ej) in equation [XTOl 
for JiE,E',ef ), 



J{E,E',ef)^{E'-ef)I{E,E',ef)- 



leading to, 



JiE,E\ef) ^ i^l^jiln 



[ef-EY 



(A22) 



r'2 (x'2 + 1) 



r2 (.t2 + 1) 



(A23) 



{E'-EY [2 
-I- j (tan^^ a; — tan^^ a;') } 
1 

^ E-E'' 



APPENDIX B: HALL EFFECT 



In this appendix, we derive equations 110611109] and [T 
in section fVIIB I Equation 11031 is 



T/(2)(±) _ 



E 



k'P±|gc/|p/)(p/|gc/|kP± 

ek-ep~{Ei-E±)+i5 



(Bl) 



Expanding for I = S using equation [ 



^\2 jsin^ 
^ p 



p)^, T i sin (fc' - p)y, sin [k - p)^, =F i sin (fc - p)^ 



ep-(i;s-i;±) + i5 



(B2) 
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Expanding the numerator inside the summation into products sinfc^, cospx' etc, 

Num = [(sinfc^/ cosp^' — cosfc^/ sinpx') T i (sinfc^/ cospyi — cosk'y, sinpj,/)] 
X [{siukx' cospx' — coskx' svapx') T i (sinfe^/ cospyi — coaky> sinpyi)]* 

I 



(B3) 



This can be further expanded into terms with products 
of the form cos^Pa;', cosp^;' cospy', cospx' sinpyi, etc. 

The denominator in equation IB2I only depends on p 
through the energy and has D4^h crystal symmetry. 
Thus, for each p, the sum includes terms with px' — > 
—Px' and Pyi — *■ —Py'- The only non-zero terms in the 
numerator have products cos^p^;', cos^py', cospx' cos Pyi, 
sin^Px', and sin^pj,'. 

The real part of the numerator is, 

Re(Num) — [ (sin /c^, sin fc^;' cos^ Pa:' 
+ cos k'^, cos kx' sin^ px') 
+ (sin k'y, sin kyi cos^ Pyi 



cos k'y, COS kyi siv? Py')] 



(B4) 



Using cos^ Px' = l/2(l + cos2p3:/) and sin^pj^' = 1/2(1" 
cos 2px> ) with similar equations for pyi , 



Re(Num) = - [cos(fc' — k)x' + cos(fc' — k)yi] 

~\. [cos(fc' + k)x' + cos(fc' + k)yi] cos2pa;', (B5) 



where we have used cos 2pj,' = cos 2pxi due to symmetry 
in the summation. The definitions of the primed coordi- 
nates in equations [551 and [Ml lead to, 

Re(Num) = cos - (fc^ — kx) cos - {k'y — /c^) 
— cos 2pxi cos i (fc^ -I- kx) cos ^ {k'y + ky) . (B6) 

A similar expansion for the imaginary part leads to, 

Im(Num) — (=p) (sin/cj;' sin/c^/ 

— sin kyi sin fc^, ) cos Px> cos Py' 

= Au'u cos Px' cos Py', (B7) 

with Ak'k defined in equation 11081 Collecting the real 
and imaginary terms. 



(2)(±) 



k'k 



X I Fo{ujo) COS - {k'^ - kx) cos - {k'y - ky) 



F2(wo) cos i (fc^ + kx) cos i {k'y + ky) 



(^^o)A'k| , 



(B8) 



where lo^) — ek ~ {Eg — E^) and the F functions are 
defined as. 



Fo{u) 



1 



Fi{lu) 



" / p - • 

1 \ cos Pa;' COS Py 



i6' 



Am 



E 



i5 



COS 2px' 

uj — Cp + iS 



Expanding the cosines in IBIOI and IBlll 



cos Pa:' cos Py' - 
COs2p;r' = COs(pa; + Py) 



1 / 

- (cos Pa; -I- COSPj,) . 



COS Px COS Py 

sin Pa; sinpj^. 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



The second term in equation IB13I averages to zero when 
summed in F2. Fi and F2 become 



Fi{uj) 



\ (cos Pa; -I- COSPy) 



i5 



^'<") - (^) T. 



COS Px COS Py 

UJ — ep + i5 



(B14) 



(B15) 



(2)(±)/ 
c'k \-^) 

with Es — > Ed and =f(z) ±{i). 

From equation 11051 the real part of V-^^^\s) is nec- 
essary for the lowest order skew-scattering. In addition, 
only terms antisymmetric under interchange of k' and k 
contribute to p^'xy'^- The last term in equation lBSl is anti- 
symmetric and only the imaginary part of Fi contributes 
to the skew-scattering. 

Defining Fi{uj) = (-l/7r)ImFi(t^), 



The equation for V^,-^ {D) is the same as equation 



Fi{u) 



leads to equations 11061 and 11071 

Equation 11231 for the skew-scattering is obtained by 
solving the linearized Boltzmann transport equation for 
skew-scattering 



9ife 
dt 



field 



dfk 
dt 



= 0, 



(B17) 
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'dfk 
dt 



field 



dfo 
dek 



i-e) (vk • E) 



(4)t^<v.AB,, (Bl«, 



^(i)(sfce») (B24) 



The current density per spin is, 



^ir) = + y^A-Wskew {k ^ k') gk 



+ Wskew {k' k) gk'] , 



(B19) 



(B25) 



where e > and fk = /o + 5fe with /o the Fermi-Dirac 
function at energy e^. l/rfc is the ordinary scattering 
rate. This is given by equation [85] in our model. From 
equations I115MTT71 J2k' '^skew(k k') = Q making the 
second term on the right hand side of equation IB19I van- 
ish. 

Substituting, 



the transport equation becomes, 



Tk ^ 



(k k')^k' = 0. 



(B20) 



(B21) 



where fi is the total volume. The 0{EB) contribution 
to the conductivity per spin (7^^'^"' from skew-scattering 

satisfies J^'^^'"'"^^ =af/^E^, 



^skew 



kk' 



dfo 
dek 



TkTk' 



X VkyVk'xWskew{k ^ k'). (B26) 
Using Uyx — —CTxy from crystal symmetry and p^y = 
Oxylo'^ for the Hall resistivity, 



„skew 2 
Pxy ^ 



kk' 



dek 

X VkyVk' skew{k ^ k ) 



(B27) 



<^k can be expanded in the series, 



^{l){skew) 
k 



(B22) 



with the additional factor of 2 on the right hand side to 
include spin. The conductivity a is, 



(l)(ord) 
fe 



is the 



where ~ 0{E) and - 0{EB). $ 
ordinary band contribution to the Hall effect. We do 
not write down the expression for this term. The second 



{l){sk&w) 



is the additional contribution 



0{EB) term, 

arising from the skew-scattering term, Wakew ■ Substitut- 
ing Ell] into [B2T] leads to. 



(0) 



-e (rfcVk • E) 



(B23) 



(-) 



m/eff 



(B28) 



where {n/m)ef / is the effective n/m and we have taken Tk 
to be independent of k and equal to t. Since py^ = —pxy 
by symmetry, we symmetrize the expression for pxy by 
setting Pxy — {l/2){pxy — Pyx), leading to equation 11231 
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